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Preface 


This is a book about prime numbers. Let 7(x) denote the number of 
primes less than or equal to x, and let Li(a) denote the logarithmic 
integral of z. According to the prime number theorem, we have 


a(x) ~ Li(a). 
This is an asymptotic statement: it means that 


1 (ax) 


> 1 5 + oO. 
hits) as £00 


If you look at any table of primes, you will find that z(x) < Li(«). 
The numerical evidence for this is quite strong. In fact, we know that 


m(a) < Li(a) 


for 2< a < 101%; see Theorem 2(1.10)]. 


Nevertheless, it remains a fact, proven by Littlewood in 1914 {Lij, 
that the sign of m(a) — Li(z) oscillates infinitely often! 


‘ 


It seems appropriate to describe this phenomenon as a “race” 
between 7(z) and Li(az). Sometimes Li(a) is ahead; sometimes 7 (x) 
is ahead. According to the numerical evidence, Li(x) starts in front. 
So, when is the first time that (a) gets ahead? That is to say, what 
is the smallest natural number for which the prime number theorem 
undercounts the number of primes? 


1X 


x Preface 


The Skewes number is, by definition, the least natural number 
with this property. The Skewes number will be denoted =. So we 
have 


(©) > Li(S) 


and = is the least number with this property. 


What is = ? No one knows. The best one can do is to provide an 
upper bound for =. The best upper bound for =, at the time of this 
writing, is a number with 317 digits! 

I should mention that there is a broad class of prime number 
races, which are pursued in the recent article by Lichtman, Mar- 
tin, and Pomerance [EMP]. Nevertheless, the race between 1(x) and 
Li(x) holds its position as the most important of these races, hence 
the title of this book. 


The mathematical background for the Skewes number is quite 
extensive, and the aim of this book is to provide this background. 
In particular, we provide proofs of the explicit formulas of analytic 
number theory, a proof of the oscillation theorem of Littlewood, a 
refinement of Lehman’s theorem, and a concise proof of the prime 
number theorem (due to Newman-Zagier). 


The main aspects of the mathematical background are 


e the Riemann explicit formula, an equation between the 
prime counting function on one side, the zeta zeros on the 
other side, 


e Lehman’s theorem, an integrated version of the Riemann 
formula, which is amenable to numerical calculation, 


e numerical analysis using a large number of zeta zeros. 


For my part, the endless fascination of the Skewes number lies in 
the interaction between the deep conceptual formulas of Riemann and 
Lehman and the subsequent numerical analysis using, for example, 22 
million zeta zeros. 

The Skewes number is named after the South African mathemati- 


cian Stanley Skewes, who was a research student of John Edensor 
Littlewood. 


Preface al 


T have written this book with an undergraduate audience in mind. 
If some of the proofs seem hard, then skip them at a first reading. To 
paraphrase the author of Ecclesiastes: there is a time to read proofs 
and a time to refrain from reading proofs. To leaven the experience 
of reading the proofs, I have provided some historical context and 
background. This includes a section on Kontsevich-Zagier periods 
and some scholarly historical material. 


Overall, this short book will serve as an invitation, a gentle in- 
troduction, to the explicit formulas of analytic number theory, using 
the Skewes number as thread and motivation. 


The book is an expansion of lectures which have proven suitable 
for third-year undergraduates in the UK. 


I take a fairly direct pathway to the explicit formulas, and in this 
context it’s a pleasure to acknowledge my debt to Jeffrey Stopple’s 


book [St]. 


Chapter 8 comprises a set of exercises which have been found 
suitable for an undergraduate course. Hints, and sometimes complete 
solutions, are also provided. 


Throughout this text, the emphasis is on the Riemann explicit 
formula. My experience is that undergraduates respond well to the 
splendor of the Riemann explicit formula. 


The Riemann explicit formula is an equation of the following 
form: 


a counting function depending only on the prime numbers 


=a formula depending only on the zeta zeros. 


Note that, in the explicit formula, the left-hand side depends only 
on the prime numbers and the right-hand side depends only on the 
zeta zeros: this is sometimes referred to as the reciprocity between the 
primes and the zeta zeros. 


This formula contains oscillatory terms corresponding to the com- 
plex zeros of the Riemann zeta function. In the book by Mazur and 
Stein [MS| p. 103], the set of the imaginary parts of the complex 
zeros of ¢(s) is called the Riemann spectrum of the primes. If each 
number in the Riemann spectrum is interpreted as a frequency, then 


xii Preface 


the Riemann spectrum is a sound spectrum, i.e., a chord in the mu- 
sical sense, with a definite pitch (the fundamental frequency) and a 
certain timbre (the overtone series). At this point in the book, there 
are some original touches: the chord of the Riemann spectrum gets 
louder as time goes by! There is a sense in which one should be able 
to hear the Riemann hypothesis. 


I have been especially influenced by Prime obsession, the excellent 
popular book by John Derbyshire [Derb], by the detailed analysis in 
Edwards of the Riemann memoir of 1859, and by Jeffrey Stop- 
ple’s book [St]. We have, above all, been influenced by the Riemann 
memoir [R]. My historical sources are 


— the articles by R. Ayoub [A], J. C. Lagarias [LE], T. J. Osler 
[Osler], V. S. Varadarajan [V], 

— the chapter on Euler in André Weil’s book [We], Euler’s book 
Introduction to analysis of the infinite [E], Riemann’s original memoir 
RI, 

— the Euler archive [Euler]. 


I first met the Skewes number in Littlewood’s Miscellany pp. 
110-112], which I read as a schoolboy. Somehow, the Skewes number 
lingered at the back of my mind, and, when in 2006-2007, I had the 
opportunity to work on it with my MSc student Kuok Fai Chao, I 
took the opportunity. This led to the article [CP], the first of recent 
papers to re-examine Lehman’s proof. 


It is a pleasure to thank Nick Gresham for help, at various times, 
with computational aspects. I wish to thank Christine Lee for many 
insightful discussions. I also wish to thank the third-year undergrad- 
uates at Manchester University who attended my course, for their 
stimulating and varied responses to my lectures. Their responses def- 
initely affected the shape of the lectures, hence the shape of this book. 


Finally, it’s a pleasure to thank Eriko Hironaka, AMS Editor, for 
her patient support and encouragement over a period of several years. 


Roger Plymen 


Chapter 1 


The Riemann zeta 
function 


1.1. Introduction 


Let z(a) denote the number of primes less than or equal to x. We 
have the logarithmic integral 
"dt 
Li(x) := —. 
9 logt 


This integral ane careful discussion. Owing to the singularity 
of the function t legit ; at t = 1, we have to define this integral via 
the Cauchy principal value taken at 1: 


1l—-e 
Li(z) := lim (/ +f a 
e0+ Lte logt’ 


The reader will sometimes encounter the function 
” dt 

li(a) :-= —. 

2 logt 


It is immediate that 


Li(x) = Li(2) + li(z) 
for « > 2. The notation for Li(x) and li(a) is not universal: the 
current notation in the Wikipedia page on the logarithmic integral 


1 
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function is the reverse of the present notation. Here, we are following 
the notation in [MV] p. 189]. 


According to the prime number theorem, we have 
a(x) ~ Li(a). 


This is an asymptotic statement: it means that 


> 1 as «© Oo. 


It is worth noting, in passing, that Li(2) + 1.045163..., so that 
m(2) < Li(2), albeit by a narrow margin. Indeed, there is a mass of 
numerical evidence which suggests that a(x) < Li(a). We have the 
following remarkable result of J. Buthe: 

n(x) < Li(a) 


for 2< x < 10'%; see Theorem 2(1.10)]. 


Nevertheless, it remains a fact, proven by Littlewood in 1914 [Lij, 
that the sign of m(a) — Li(x) oscillates infinitely often! 


Consider the set 
{xz EN: x(x) > Li(x)}. 


Littlewood’s theorem implies that this set has a least element. 
The least element in this set of natural numbers will be called the 
Skewes number. We will denote the Skewes number by =: 


E:=inf{e ©N: r(2) > Li(x)}. 


We can be sure that 10!° is a lower bound for the Skewes number 


(0 


10% < =. 


What is =? No one knows. The best one can do, at the time of 
writing, is to find an upper bound for =. Such an upper bound will 
be called a Skewes number. 


In due course, Skewes established his famous upper bound: 


103 
exp exp exp exp(7.705) < 19!" 
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Theorem 1.1 (Skewes {S]). For some x < exp exp exp exp(7.705) we 
have 
n(x) > Li(z). 


In this remarkable article, the proofs are based on both of the ex- 
plicit formulas: the von Mangoldt formula (see Section 4.1) is quoted 
in [S| Part I, Eq. (2)], and the Riemann formula (see Section 4.6) is 
quoted in [S] Part II, Eq. (24)]. 


Subsequent developments by Lehman and others all depend cru- 
cially on the Riemann explicit formula. 


How do we reduce the upper bound? Lehman decided to integrate 
the difference a(x) — Li(az) against a positive test function over a 
certain interval [a,b]. If the resulting integral is positive, then we can 
infer that a(x) > Li(a) for some x in the interval [a,b]. To prove 
that the Lehman integral is positive, the only known way is to enter 
millions of zeta zeros with the help of a powerful computer. This 
method has been successively refined, and the current best known 
bounds are as follows. We will write them in the natural log scale: 


19 log 10 < log= < 727.951346801. 


There are some very large natural numbers, for example, the 
estimated number of stars in the universe, or the estimated number 
of viruses in the oceans of the world, but these numbers are dwarfed 
by the upper bound exp(727.951346801), which is a colossal number 
with 317 digits. 


1.2. The Riemann zeta function 


In this section, we will meet the Riemann zeta function, named after 
Bernhard Riemann (1826-1866). The notation ¢ for the function and 
s for the variable are both due to Riemann in his memoir of 1859. 


Euler (1707-1783) had worked extensively on infinite series of the 


form 
[oe 


1 
Lo ae 
where k is a positive integer, see his famous book Introduction to 
analysis of the infinite, Chapter X [E]. In Riemann’s notation, we 
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would write 
Co 


1 
C(k) = Do 
n=1 
In fact, Euler discovered an elegant formula for ¢(k) where k is an 
even positive integer. 


Riemann’s inaugural dissertation of 1851 was called “Foundations 
for a general theory of functions of a complex variable” . With this in 
his mathematical background, it was perhaps natural for Riemann to 
interpolate between the positive integer values of k, thereby leading 
to the zeta function ¢(s) of a complex variable s. With hindsight, this 
may seem like a natural, even inevitable, step, but it is surely a big 
conceptual leap. When this is done, there are issues of convergence 
and we have 


with R(s) > 1. 

We will show in this chapter that the zeta function can be ex- 
tended to a meromorphic function in the entire plane C. We briefly 
review some of Euler’s pathbreaking work on ¢(2k). 


With acknowledgement to Euler, the zeta function ¢(s) is some- 
times called the Euler-Riemann zeta function. 


Riemann’s memoir is to be found in [R], and English translations 
are readily available; see for example [Ed]. 


1.3. The prime numbers 


Definition 1.2. An integer n > 2 is prime when the only divisors 
of n are 1 and n itself. If n > 2 is not prime, then it is said to be 
composite, or nonprime. 


We can list the first few prime numbers: 
2,3, 5, 7,11, 13,17, 19, 23, 29, 31,37.... 


Then given any positive integer x € Z we define z (x) to be the num- 
ber of prime numbers less than or equal to x. For example x (10) = 4, 
m (20) = 8, 7 (40) = 12... 
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Theorem 1.3. There exist infinitely many prime numbers. 


Proof. This is from Euclid’s Elements IX 20. It is really a one-line 
proof: let p1,p2,..-,Pn be the first n prime numbers and consider 


Q=1+ pi: p2--+Pn- 


To spell out the details: 

— EITHER @Q is prime, in which case there is a prime, namely Q, 
greater than py, 

— OR Q is nonprime, in which case each prime factor of Q must 
be greater than each p,;, for division of Q by p; admits remainder 1. 

Either way, there is a prime greater than p,. This proof does not 
generate primes, but it does place an upper bound, namely Q, on the 
next prime after p,. The first few steps yield 


1+2=3, 
14+2-3=7, 
14+2-3-5=31, 


1932527 =211, 
142-3-5-7-11 = 2311, 
142-3-5-7-11-13 = 30031. 


Of these, 3,7,31,211,2311 are prime; 30031 is composite with 
least prime factor 59. 


Definition 1.4. A Mersenne prime is a prime number of the form 
M,, = 2” —1 for some integer n. 


They are named after Marin Mersenne, a French Minim friar, 
who studied them in the early 17th century. As of December 2018, 
the largest known Mersenne prime is 


982,589,933 <2 1, 


a number with 24, 862, 048 digits. 
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Note the algebraic identity 
ola“ Dilys tetety™*) 


and set y = 2”. Then we have that 2” — 1 divides 2” — 1, and so if 
n € Z is composite, then 2” — 1 is composite. In other words, if M,, 
is prime, then n must be prime. Each Mersenne prime is therefore of 
the form 

M, = 2? —-1. 


For more information on this topic, go to GIMPS, the Great 
Internet Mersenne Prime Search. 


1.4. The Riemann zeta function 


We begin with the definition of the zeta function. This definition first 
appeared in the famous Riemann memoir of 1859. 


Definition 1.5. The Riemann zeta function is 


co 
1 
C(s) = ae 
n=1 
with s eC, R(s) > 1. 
Writing s = 0 + it, we obtain 
1 —s 
nm 
- e slosn 
os e 7 log n—itlogn 


_ e 7 legn : ett logn 


and so ; 
— =n ° (cos(tlogn) — isin (tlogn)) . 
n> 
We can now write 
“= cos (tlog n sin (tlogn 
R(c(s)) = Jo EBM gna (G(s)) =— So BH) 
n=1 n=1 


and also |1/n*| = 1/n?. 


It’s useful to study the zeta function for some real values of the 
variable s. 
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oo 1 
n=1 n> 


Theorem 1.6. The infinite series given by > has the following 


properties: 
(1) The series converges if s > 1. 
(2) The series diverges if s =1. 


Proof. Suppose first that s > 1; then observe that 


1 1 1 1 2 
| gl-s 
+ < 58 < 38 = 


and 


4s 5s 65 7s 4s 4s 4s 4s 
and so on. The geometric series 


Co 


f(s):= oa)" 


n=0 
converges if s > 1 for then we have 0 < 2!~* < 1. This gives us that 
s>l=> s = <0 
n=1 uy 


Next, we consider the case when s = 1 and we consider the nth 
harmonic number 


For the harmonic numbers we have the following elegant and im- 
portant inequalities: 


Lp tard fd 
Cee nae Wee ae, 
re ee ee ee i oe oe 
5 be FOB By BBR ~ De 
ne rae ane ee ee ite gale eae 
9' 10° 11 PAG A6e 1616. 16 2” 
and so on. In general, we have 
and so 


Hy: +00 as ko. 
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The conclusion is that the harmonic series 
1+ ! + : + : + 
2 3 +4 


is a divergent series. 


The divergence is surprisingly slow. Here’s a thought experiment 
to emphasize this point. A computer is programmed to add 1 trillion 
terms of the harmonic series )> A each second. Since the Big Bang, 
has enough time elapsed for the Nth partial sum to exceed 100 ? 


Since 
| 
ns|— n®(s)’ 
it is clear that 
Co 
1 
ne 
n=1 


is absolutely convergent in the open half-plane R(s) > 1. 


1.5. Euler and the zeta function 


Euler found formulas for ¢(s) when s is an even natural number. 
In his book Introduction to analysis of the infinite, he showed that 
pp. 139-140] 


Sd Oy 
ea 
jy ee 

~ 24 nt ~ 5 Ty 


ns 
n=1 
a) 2 aa 
¢(10) = >) 5 = yr 3™ ’ 
n=1 


1 2!691 45 
<)>) aati” 
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a ll 212 35 
14 = Wns Seta 14 
(19 =) oat 
1 243617 4, 
616) =) ae Te 
eel 216 43867 
18) = - 2 18 
cue) pater i9) 21’ 


1 218122277 45 
6(20) = dU no O11 6B ° 


n=1 

=. - BPO B54519 4, 
((22) =) om = oar gs 

n=1 


3 1 2? 1181820455 54 
m4 951 273~— 


=I 
“. 1 27476977927 96 
(26) =) eo 


For the original document, see E130|1739, §24]. Return- 
ing to his book, Euler writes p. 140): 


...we have gone far enough to see a sequence which at first seems 
quite irregular, 1,1/3, 1/3, 3/5, 5/3, 691/105, 35/1,..., but tt is of ex- 


traordinary usefulness in several places. 


Later, Euler discovered the pattern underlying these apparently 
irregular numbers and obtained his beautiful formula 


_ (—1)*+122*-1 Bo, . 
eek) = (2k)! ; 


Using 7(s) = (1 — 2!~*)¢(s) we have 
n(2k) = 3 oe Ale = ay aah. 


nk (2k)! 


n=1 
Here the numbers B;, are the Bernoulli numbers, and they are all 
rational. For the first few we have Bp = 1, By = —1/2, and 


B3 = By = By =::--=0 
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and their generating function is 


Notice that ¢(s), with s an even natural number, is a transcen- 
dental number. 


1.6. Meromorphic continuation of ¢(s) 


We say that a function f(s) with s € C is holomorphic in an open set 
U CC if the derivative df/ds exists for all s € U. We write Hol(U) 
for the set of all complex functions that are holomorphic in the open 
set U. The set Hol(U) is a vector space: for any f,g € Hol(U) and 
EC, we have 


d(f+g) _ df es dg 


ds ds ds’ 
arf) _, (af 
ds ds }~ 


The kind of convergence which is useful in complex analysis is 
rather strict: we require uniform convergence on compact sets. Uni- 
form convergence on compact sets is called locally uniform conver- 
gence. The usefulness of this kind of convergence is spelled out in the 
following important lemma. 


Lemma 1.7 (The convergence lemma). Let U C C be an open 
set, and let f, € Hol(U) be a sequence of functions. If fr —> f 
locally uniformly, then f € Hol(U) and 
dfn df 
Pani ail —_ pataliey 
ds ds 


also locally uniformly. 


as n— oo 


We continue with an immediate application of the convergence 
lemma. 


Theorem 1.8. The Riemann zeta function ¢(s) is holomorphic in 
the open half-plane R(s) > 1. 
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Proof. Let D = {s € C: R(s) > 1}. Note that for any 6 > 0, for all 
s €D such that R(s) > 1+, we have 


1 


n> 


1 
ae 


In other words, since for each n the function 1/n* is holomorphic 
everywhere and the sum converges uniformly on compact sets of D, 


¢ is holomorphic on D by the convergence lemma. 


Before discussing the meromorphic continuation of the zeta func- 
tion, we recall the four types of isolated singularities: removable sin- 
gularity, simple pole, pole of order n, essential singularity. We illus- 
trate via singularities at 0. 


(1) Removable singularity at 0. 


snz 1 gent 
z Toes ) (2n + 1)! 
ze gt 6 
Pe get oy a 


The apparent singularity at z = 0 can be removed. 
(2) Simple pole at 0. 


oe) 
y2ntl 


sin z pa 
s Y) (2n + 1)! 


1 
2 7 


z 
n= 


1 z Zz z 


The 1/z term gives rise to the simple pole at z = 0. 

(3) Pole of order n at 0. An example of a pole of order n is an 
expression of the form 1/z”. This is a pole of order n at z = 0. 
Note that a simple pole is a pole of order 1. 
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(4) Essential singularity at 0. Here is a good example. 


spk a cl 
oy Foe} * 3th 28 , 


Definition 1.9. A function f(z) of a complex variable z € C is mero- 
morphic in an open set D C C if it is holomorphic at all points except 
for isolated poles. The function f(z) is entire if it is holomorphic at 
all points z € C. 


Now consider the power series 


provided that |z| < 1. The function f(z) = 1/(1 — z) is defined 
everywhere in C — {1} even though the series }~°*_) 2” only converges 
in the open disc |z| < 1. 


Note that f(z) = 1/(1 — z) is not defined at z = 1 as there is a 
simple pole there. But by Definition (1.9), this function is meromor- 
phic as the pole is an isolated singularity. This idea is called mero- 
morphic continuation, which we now illustrate with the Riemann zeta 
function. 


Theorem 1.10. The meromorphic continuation of the Riemann zeta 
function ¢(s) into the open half-plane R(s) > 0 is given by 


(1.1) ((s) =1+ sf Ao ae. 
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Proof. For §t(s) > 1 we have 


where [¢] denotes the integer part of the real number t and {t} denotes 
its fractional part. 

Since {t} € [0,1], the last integral converges for R(s) > 0. The 
right-hand side thus defines a meromorphic continuation of ¢(s) to 
the half-plane R(s) > 0. 


The simple pole at s = 1 is clearly visible in equation (LI). 


In order to extend the domain of the Riemann zeta function to 
the entire complex plane by the method of analytic continuation, 
we need to adopt a more systematic approach as in [D] 15.5]. This 
is accomplished step by step by extending the function first to the 
domain Do = {z € C: R(z) > O}, next to the domain D_; = 
{z €C: R(z) > —1}, and so on, adding a vertical strip to the domain 
at each step. 


Observe that for each n, the function gn(z) = — meet is 
entire, with a zero at 1. Hence, for each n, h,, defined by 


os 


is also entire. But for each n, 


Me igiy ma ie 
(1.2) m= f oSo= fs. 
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Similarly, the function f,,(z) defined by 


1 
fn(z) = @+i hn(z) 


od [ dt 
~ (n +1)? 9 (n+1)? 
a tdt 
=-Zz 
a Gore Leet 


is also entire being the difference of two entire functions, as may be 
verified directly by integrating the last term by parts. 


Now, for each z, 


‘ tdt 
|fn(z)| < || f Imtoee 


2 . tdt 
< |z|- 5 Mt pR@H 


and if z is such that R(z) > 6 and |z| < M for some 6 > Oand M > 0, 
we have 


1 
lfn(z)| <M - PEE) 


so that 5> f, converges uniformly on compact sets of Do and hence 
defines a holomorphic function on Do. Note that while each f,, is 
entire, >> fp is holomorphic on Do only. 


On the other hand, we have for z € Dj, 


(=1l+e+at 
= 1 
= eens 


=1I+ > (/ (n = 7 a [ (n =) 


n=1 


(7). ae =f tae 
=1+4+ a a a —_ # 
ay (n+t)? Lf (nto 


n=1 


Using equation (1.2), we have 
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Thus, 


(1.3) ¢(z) = 


coef oat 


n=1 


Thus, for z € Dy, ¢(z) is given by the above expression. But in this 
expression the last term is a holomorphic function on Do while the 
function 1 + 1/(z — 1) is meromorphic on Do with a simple pole at 
1. Hence the function defined by the right-hand side of is a 
meromorphic function on Do with a simple pole at 1, and it coincides 
with the definition of the Riemann zeta function on D,;. Hence we 
are justified if we define the zeta function on Do by (1.3) above. This 
extends the zeta function from D; to Do. 


We proceed to continue the function to the domain 
-1={zeEC: R(z) > -I}. 


The explicit formula is 


1 z ae dt 
(a) =14+ + - Fie ti)-] oe wes 


which is now holomorphic everywhere on D_, except at 1, where it 


has a simple pole. 


The next explicit formula is 


z—-1 2! ! 
(n+ t)2+3 


which is now holomorphic on D_2 = {z € C: R(z) > —2}, except of 
course at s = l. 


The technique to go from one explicit formula to the next is little 
more than integration by parts. 


In this fashion, we define the Riemann zeta function over the 
entire complex plane. It is a meromorphic function, with a single, 
simple pole at z = 1. The residue of the zeta function at z = 1 is 1. 
Note that, at each step in the proof, the exponent of (n + t) in the 
denominator of the integrand increases by 1. 
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Definition 1.11. The number p € C is a zeta zero if ¢(p) = 0. 


The Riemann hypothesis. All the zeta zeros in the open half-plane 
R(s) > O lie on the line R(s) = 1/2. 


The Riemann hypothesis is at present unproven. It is one of 
the most famous problems in analytic number theory, which holds 
the key to understanding the distribution of the prime numbers. 
The Riemann hypothesis is one of the Millenium Problems; see www. 
claymath.org/millenium-problens. 


Chapter 2 


The Euler product 


2.1. The zeta function and the Euler product 


In this chapter, we meet the Euler product 


with R(s) > 1. This is an infinite product taken over all the prime 
numbers p = 2,3,5,7,.... We shall see that the Euler product is equal 
to ¢(s). This depends ultimately on the uniqueness of factorization 
for the natural numbers. 

The Euler product can be traced back to Euler’s seminal book 
Introduction to Analysis of the Infinite Chapter XV]. In that work, 
Euler takes s to be a positive integer. The product symbol [] was 
not available to Euler, and he wrote 


(Gael ea 


The Euler product is not immediately useful; however the log- 
arithm of the Euler product is crucial. The logarithm of the Euler 
infinite product is an infinite series, the sum taken over all the prime 
numbers p. This infinite series is exploited later in this chapter to 
prove the important result that 


R(s) = 1 = > C(s) £0. 


18 2. The Euler product 


The logarithm of the Euler product is further exploited, later in 
this book (see Section 4.7) in proving the Riemann explicit formula. 


For an informal proof of the equality of the Euler product with 
the zeta function, we can proceed as follows: 


Pp Pp 


al 
a 


I 


=1 


¢(s) 


However, in order to take care of convergence issues, we may 
proceed as follows. We can think of the following theorem as the 
first theorem in analytic number theory. Note that Euler already had 
this theorem (with a forgivable typo) when s is a positive integer, see 


§275]. 


Theorem 2.1. The Riemann zeta function ¢(s) is equal to the Euler 


product; that is, 
¢(s)=]](l-p) 
Pp 


with R(s) > 1. 


Proof. Recall that for R(s) > 1, 


1 1 1 
SS ed EE Sadi 
Oe ne ce eee 
Therefore we have 
1 1 1 1 1 
i =] sega fe cseace 
( =) st) 3s T 5s tas TF Gays 


as all the even terms cancel each other. Then 


(5) (-B)09-(-S) Gem eeRen) 


iL 1 1 
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This time all terms having 3 as factor in the denominator cancel each 
other; this is an example of the sieve of Eratosthenes. Proceeding in 
this way, if p1,p2,...,px are the first k primes, 


: i 1 
TI (1-55) s- = -- Wie 


with summation over all natural numbers m which do not have p1, 
p2,--+;Dk as a factor. 


Let gz be the least such number, and let R(s) = 1+6 with 6 > 0. 
Then we have 


: 1 1 
fi(i-8)enaf eS gta 
j=l P3 m>qn 


The term on the right-hand side is the tail of the infinite series for 
¢(1+ 6). Now let k — oo. We get 


H(:- 3) <= 


and the convergence is locally uniform in #(s) > 1. This gives Euler’s 
theorem for R(s) > 1 


where the product on the right is over all the primes. 


This proof provides us with a bonus, in the following way. 


Theorem 2.2. The Riemann zeta function ¢(s) 4 0 in the open 
half-plane R(s) > 1. 


Proof. Given 0 < « < 1, there exists k such that 


1 
(1-5) ce -1 <€, 


j=l j 


This inequality forces ¢ to be nonzero in §(s) > 1. 
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2.2. The logarithmic derivative of ¢(s) 


We now consider the logarithmic derivative of the Riemann zeta func- 
tion. First, we need to consider what the derivative of an infinite 
product is. Suppose that f is a function that can be written as 


i= [a 
i=l 


So what can we say about the derivative df/dxz? To get an idea 
how this works, let’s first consider the derivative of finite products of 
functions. By the familiar Leibniz rule for products, we have 


d(gig2) _ dg dgz 


de dae de 
and 
d(gigegs) _ dgi dg2 | dgs 
de 9298 a + 9193 Ae 9192 re 
So we have 
f = 919293 


1 df = yl dgy , 1 dg2 ‘ 1 dg3 
=7lz)-a(2) | ala) | mes. 


Similarly, for a finite product we get that 


1 (df “1 (dg 
=>>—— | =—s = aoe atts : 
jae) “2a ae) 
This is a good pointer as to what happens for infinite products. We 


now move on to infinite products. 


Theorem 2.3. Let U be an open set in C. Suppose that 


f= [9 
w=1 


with g; © Hol(U), f £0, g; 40 in U for alli > 1, and suppose that 
the convergence is locally uniform in U. Then we have 


1(df\ Sl (dg 
(ge) Deg (ee) 


i=1 
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Proof. Let fn = 4 g; be the product of the the first n functions. 


We have that 
1 (dfn\ wr (dg 
aca lree 


By the convergence lemma, we know that f € Hol(U) and that 
df, /dz —> df/dz asn —> cw. Let C CU bea circle. Since C is 
a compact subset of U and |f(z)| > 0 on C, there exists 6 > 0 such 
that |f(z)| > 6 > 0 for all z € C. Let N be such that |f, (z)| > 6/2 
for all zg € C and for n > N. Then, for z€ C andn> WN, 


_ lin (2) - F 
If @) FO 


< Bln (2) -F (2). 


1 1 
fn(2) fF () 


So 1/fn — 1/f uniformly on the circle C. Hence, by the maximum 
modulus theorem p. 144] we have that 1/f;, —> 1/f uniformly 
on every compact disc D C U. 

At this stage of the argument, we know that 1/f, converges 
locally uniformly to 1/f and df,/dz converges locally uniformly to 
df /dz. Combining these two results, we infer that the sequence f// fr 
converges locally uniformly to f’/f. 


Hence 


as required. 
If we have a holomorphic function F(z) in a domain U, then we 
can differentiate after taking the logarithm: 


(log F(@)) = FI. 


This expression is called the logarithmic derivative of F’. 
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The reward for our recent rather technical excursion is the fol- 
lowing result on the logarithmic derivative of the zeta function. 


Lemma 2.4. For the logarithmic derivative of the zeta function, we 
have, for R(s) > 1, 


aa (#) = a5 (&(He--9")) 


Pp 


Bee 


(l—p-*)* 
eS (1—p7*) ? p~*logp 
ie , 
3 (tp) 

2 -S- (1 —p~*) ‘pS logp 
Pp 

=e p *logp 
eae) 

= ->- log p 
Pp ” 7 1 


We now introduce the von Mangoldt function. 


Definition 2.5. We define the von Mangoldt function A(n) as fol- 
lows: 


A(n)=logp if n isa power of p 


=0 otherwise. 


Theorem 2.6. The logarithmic derivative of the Riemann zeta func- 
tion in the open half-plane R(s) > 1 is given by 


C(s)_ An) 
Ray 2 ae 
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Proof. Suppose that s € C and consider the zeta function ¢(s) in 
the open half-plane #R(s) > 1. Now by Lemma[2.4] we have that 


Ci(s) _ p* 
Ge) ~ 2 ee a 
=) 7 -logp op * 
p k=1 


= ~log2S2-** ~1og3 3 ~log5 55 = 
k=1 k=1 k=1 


Ss Al(n 
a> och 


This completes the proof of the theorem. 


We know that the series }7°°_, 4 is divergent. Suppose that we 
delete all composite numbers n, leaving the series 2 . where p runs 
through all the prime numbers. Is this series still divergent? 


Theorem 2.7. The series a : is divergent; that is, 


Proof. Inspired by [TMF] p. 13], write n = gm? where q denotes a 
square-free integer. We have 


~ 1 <1 ~ 1 
a do. tae 
n<ax qSt © m<y/u/q 

1 1 
< — = 
er 


<2.0-. 


qgu 
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We also have 


qx pox 
1 
< exp oS = 
psa ss 
and 
1 n+1 dt 
= a 
ee 
n<ax n<a 
> log x. 


We neatly combine all three inequalities: 


1 
S- — > log log x — log 2 


pSu 
for all « > 2. This implies that a 1/p=o. 


So we do indeed have that the series ae 1/p = o0 is divergent, 
but the divergence is very slow. Euler would probably have written 


1 
S- — = loglog oo. 
Pp 


Theorem 2.8. Suppose thats =a+it € C. The Riemann zeta 
function ¢ (s) £0 on the line o = 1. 


Proof. Suppose that s = o + it and consider the zeta function ¢(s) 
in the open half-plane 0 > 1. So we have that 


log¢(s) = S “log (1— p*) 


s 


and so # (log ¢(s)) is given by 
m7 me s(t l — =p Ne: S tl ae . 
S s P cos (tm log p) 2 y 7? cos (t log p™) 


p m=1 
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We now use a simple trigonometric identity: 


3+ 4cosy + cos2y = 3+ 4cosy + 2cos? y— 1 


=2+4cos y+ 2cos* y 
= 2(1+ cosy)? 
> 0. 


Set y = tlogp™ and we immediately have 
3 + 4cos (tlog p™) + cos (2t log p™) > 0. 


Now we sum all these nonnegative terms: 


S- 2 = (3 + 4cos(t log p™) + cos(2t log p™)) > 0 


p m=1 


which is equal to 
3 - log ¢(o) + 48 (log ¢ (o + it)) + R (log ¢ (o + 2it)) > 0. 
This statement is equivalent to the following inequality: 
# (a) (C4 (o + it)| |¢ (o + 2it)| > 1. 
This depends on the fact that 
RK(log z) = log |z| 
which follows from the definition of log z: 
log z = log |z| + ¢arg z. 


At this point, we can be sure that 


(0. i@=1)3 


1 
co) |e Saas Nic IG (o + 2i8)| > —. 
= 
Now suppose that ¢ (1+ it) = 0, with t 4 0. Recall that ¢(s) has a 
simple pole at s = 1 and so (0 — 1)¢(a) > 1 as o —> 1 from the 
right. 
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Now let o —> 1 from the right. We have 
(@-1 oe) —*1, 
IC (o + 2it)| — [C(1 + 244)], 
¢(o + it) —¢ (1+ it) 
a—1 


> O(1 + it), 


and so, as 0 —+ 1 from the right, the left-hand side of (2.1) tends 
to a finite limit. This contradicts the fact that the right-hand side of 
(2.1) tends to oo. Hence ¢ (s) 4 0 on the line o = 1 as required. 


ee |= 


We now make an elementary but slightly subtle observation. The 
critical strip is defined as {s € C : 0 < #(s) < 1}. Suppose that 
we take a rectangle D, of finite area, in the critical strip, with its 
right-hand edge on the line R(s) = 1. There will only be a finite 
number of zeta zeros within this rectangle. By narrowing the width 
of this rectangle to avoid all of these zeros, we can construct a new 
rectangle E, also of finite area, in the critical strip, its right-hand edge 
on #(s) = 1, which is completely free of zeta zeros. 


In this precise sense, Theorem [2.8] allows us to extend the zero- 
free region of the zeta function into the critical strip. 


Chapter 3 


The functional equation 


We revisit the zeta function in the open half-plane §t(s) > 1. In this 
domain, the zeta function has an evident reflection symmetry: 


cose 


co 
n=1 

lee) 
n=1 


=| = 


= ¢(5) 
where the overline denotes complex conjugation. 


The meromorphic continuation of the Riemann zeta function ¢(s) 
into the open half-plane #(s) > 0 is given by 


¢(s) =1+ : sf ae 


s—l1 


see Theorem [1.10} Taking complex conjugates, we obtain 


ja1t es [Shae 


= ¢(8). 


The critical strip is defined to be the open set 


{s:0< R(s) < 1} 
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and the critical line is defined to be the line 
{s: R(s) = 1/2}. 


So this reflection symmetry continues into the critical strip. We 
have 


¢(s) =0 = > ¢(3) =0. 


In the critical strip, the zeta zeros occur in complex conjugate pairs. 
In particular, if p = 1/2 + 7y is a zeta zero on the critical line, then 
p = 1/2 — 77 is also a zeta zero on the critical line. 

The zeta function has, in addition to this symmetry, a hidden 
symmetry, which only becomes evident after an overhaul of the zeta 
function. To carry out this overhaul, we have to multiply ¢(s) by 


Ss 


sels —1r (5) T 


tol 


where I'(s) is the gamma function. The resulting function will be 
denoted €(s). This new function has the following symmetry: 


(3.1) €(s) = €(1— 8). 


The map s +> 1 — s is reflection in the horizontal line S(s) = 0 
followed by reflection in the critical line R(s) = 1/2. 


The equation B.1) is the functional equation of Riemann, and 
appears in his memoir [R] of 1859. The functional equation was 
anticipated by Euler. Euler’s paper, written in 1749, appeared in 
1768 in the Memoirs of the Berlin Academy; it was never reprinted. 
For a hundred years after their discovery, Euler’s functional equations 
were utterly forgotten; see p. 276]. 


In this chapter, we will do many things. We will define the gamma 
function I'(s) and establish its basic properties. We will then establish 
the Euler-Riemann functional equation, using it to locate the trivial 
zeros and to compute a few zeta values. We will then give an account 
of Euler’s great work of 1749. 
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We will encounter the Euler constant, which reappears in several 
places: 


e in the definition of the gamma function; 

e as the constant term in the Laurent expansion of ¢(s) around 
s=1; 

e in the linear term of the Taylor expansion of I'(s) around s = 1. 


We make an excursion into the realm of the Kontsevich-Zagier 
periods, with the aim of providing a context for the Euler constant. 


3.1. The gamma function 


Definition 3.1. If m is a natural number, then n! is defined as 
nl =1-2-3---n. 
To this we add the convention that 0! = 1. 
Recall from Chapter 1 that the Riemann zeta function started 


life as ¢(k) where & € N. In other words, it was originally defined by 
Euler as 


We then extended the definition to ¢(s) where s € C. To define the 
gamma function I'(s), what we actually do is extend the definition of 
the factorial function n! over the complex numbers s € C. 


First of all, we define 


Tn =| et" dt 
0 


with n € N. Integrating by parts, we have 


[oe) 
i [-e*e"] > +f nett” dt 
0 


Since 


we can agree that 
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This equation is a stepping stone towards a definition of the gamma 
function. 


Definition 3.2. Given z € C, the gamma function I'(z) is defined 
by the integral 


T(z) = i. a teed. 
0 
Note that [(n +1) =n! 


The integral 


CO 
| et? 1dt 
0 


is convergent if #(z) > 1 and defines a holomorphic function in that 
domain. Let 


D, = {z: R(z) > n}. 
With z € D, we have [ (z+ 1) = 2T'(z) and so 


jones 


z 


This leads us to define the gamma function in Do as follows: 


T(z) := Biers), z€Do. 


z& 


Iterating this, we define 


T(z) := Bea: z€D-i, 


z 


and so on through all the domains D_; C D_,z C D_3C.... 


In this iterative procedure, do we encounter any singularities? 
We have 


r= f eet ae 
0 


1 fore) 
=f was [ et? dt 
0 1 


=T; (2) +L (2) 
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where 


1 
T; (z) -| As ames 
0 


and 
T2 (z) -| ae alg 2 
1 


Now [2 (z) is an entire function; in other words it is holomorphic 
in the whole of the complex plane. So the singularities of the gamma 
function are the singularities of Ty (z): 


1 
T; (2) =} et? ldt 
0 


ace 1) 
= ae ee 
> ” pztn 
Belntere) (z+n) a 
= 5 n! ( (z 2 n) 
So we have simple poles at 
z=0,—-1,-2,-3,.... 
Theorem 3.3. Each term in the above expression reveals a simple 


pole of the gamma function T (z) at the points 0, -1, —2,...,-— 
the residue at z = —n is (—1)” /nl. 


ye and 


Restricting to a real variable x, we still have 


T(a +1) 


I(x) = 


and it is instructive to draw the graph of y =I (a) for all x ER. 
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For example, let ¢ be small and positive. Then we have 


re) = 049 2 


which is large and positive. At the same time, 


r( pee Ee 1 


—€ —€ 
which is large and negative. The gamma function y = I(x) has a 
severe discontinuity at = 0, namely a jump discontuity with infinite 
jump. 


In general, and for the same reason, the gamma function y = I'(z) 
has jump discontinuities with infinite jump at each negative integer 
ee ee: Se 

In addition, the vertical lines {s : R(s) = n} with n = 0, —-1, —2, 
—3,... are asymptotes of the curve y = T(z). 

The simple poles of T(z), which occur at z = 0,—1,—2,—3,..., 
correspond to the jump discontinuities of the curve y = (a), which 
occur at x = 0,—-1,—-2,—-3,.... 

The point of view of Emil Artin is “that the gamma function I(x) 
can be thought of as one of the elementary functions, and that all of 


its basic properties can be established using elementary methods of 
the calculus.” See his monograph [Art]. 


3.1.1. An alternative approach. This approach is due to Desh- 
pande p. 225]. Consider the function ¢ defined by the infinite 
product 


co 


o(z) = [[G + 2/3) exp(-z/4). 
j=l 
This infinite product satisfies the hypothesis of the Weierstrass factor- 
ization theorem Prop. 15.14] and so defines an entire function 


wtih simple zeros at —1,—2,—3,.... The function w(z) = z¢@(z) is 
also entire with simple zeros at 0,—1,—2,.... Let ¢, denote the finite 
product 


k 
ox(z) = [ [+ 2/9) exp(—2/4). 


gat 
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It is an entire function for each k = 1,2,... and lim¢d, = ¢. Routine 
computation shows that 


#2) _ Sr A/n) — (1/n)(. + z/n) 
(z) = 1+ (z/n) 
E 1 1 
- 2d (— = | 


Similarly, 
¢%(z-) _< 1 1 
ReaD cree 
so that 
aes #(z-1)_ 2)_1_ 1 


ou(z—1) bez) 2 k+z 
We use this identity in the next result. 


Lemma 3.4. Let ¢ be the function defined above. Then there is a 
constant y such that for any z, 


o(z — 1) = exp(y) - 26(z). 


Proof. The function y¥(z) = z@(z) is entire with simple zeros at 
0, -1, —2, —3,.... So is the function w2(z) = ¢(z — 1) with the same 
simple zeros. Hence by the Weierstrass factorization theorem there is 
an entire function y such that $(z — 1) = exp(7(z)) - z@(z). We now 
show that 7 is a constant function. For this, it suffices to show that 
7 =0. 
Differentiating both sides of ¢(z — 1) = exp(7(z)) - zdé(z), we get 
'(z — 1) = exp(y(z))(O(z) + 26'(2) + ¥(2)20(z)), 
eG face) Nem aes ck ae 
z-1) z° Az) | 


From this we get 


(3.3) y(2) = 
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On the other hand, we have by (8.2) 


Pu(Z—-1) lz) _ 1 1 


de(z—1) ge(z) 2 k+z 


Now differentiation is a continuous operator on entire functions. 
bea fe : 
Hence lim ¢), = ¢' since lim ¢, = @. Hence 


w(Z— 1) Pelz)\ _ isi 
ae (SH = 2 eae € as) 
1 
_#2-1)_ #@) 
o(z-1) oz) 


From (3.3) we obtain immediately that 


7(z) =0 


so that y is a constant function. 


The constant y so obtained is known as the Euler constant. This 
is its first appearance in this book. To evaluate it, note that ¢(0) = 1 
by the definition of ¢. 


We have 
(1) = [ [G+ 1/s) exp(-1/J) 
j=l 
= limc, 
where 


n 


en = [] (0+ 1/3) exp(-1/9) 


n n 


= exp - 0/5) [[G+0/. 


3.1. The gamma function 35 


This gives 


log en = — S°(1/J) +S°( (log(j + 1) — log j), 
j=1 


1 1 
— log cp = (+5 ates -) logn — log((n + 1)/n), 


: ; 1 1 
~ Jim logen = Tim. (145 45-+2-toen). 


Letting z = 1 in the equation $(z — 1) = exp(7) - z@(z), we have 


(0) = exp(y) - ¢(1), 
log (0) = 7 + log 9(1), 
0=7+ log(limc,), 


y =-—limlogcn, 


so that 


= li beet ees 
peer Rag geodetic ogn}. 


Recall that @¢ is defined by 
o(z) = IIa + 2/j) exp(—z/j). 


It has simple zeros at —1,—2,—3,... and it satisfies the relation 
(2-1) = €724(2). 
The gamma function is now defined by 


exp(—72) 
T(z) = ——_— 
= 6) 
for z # 0,—1,—2,... where y is the Euler constant. It is clearly a 


meromorphic function with simple poles at 0,—1,—2 
ZeYOS. 


,--. It has no 


From the recurrence relation of ¢, we have 


(z+ Io(z +1) =e 762). 
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Substituting this in 
e7 (2+) 


~ (g+1)d24+1) 


gives 


which is the recurrence relation for I. 


The gamma function T(z) is a meromorphic function in the com- 
plex plane C with simple poles at 0, -—1, -2,—3,.... The gamma func- 
tion is zero-free; i.e., it never takes the value 0. We have T(n+1) = 
nt. 


3.2. The functional equation 


Recall that the zeta function started life, in the work of Euler, as a 
function of natural numbers; that is, ¢ (k) = 07°, +. This definition 
was then extended to the absolutely convergent series 


1 
= — R 1 
Oe De (s) >, 
which was then extended to a meromorphic function in the whole of 


the complex plane, with a simple pole at s = 1. 


Then at the beginning of this section we saw that the gamma 
function started life as the factorial function n!. This can then be 
extended as a meromorphic function into the whole of the complex 
plane, with simple poles at the negative integers 0,—1,—2,—3,.... 
These two functions are linked by a single equation called the func- 
tional equation, and from this we can uncover the symmetries of the 
zeta Zeros. 


The following theorem is one of the most important theorems in 
mathematics. 
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Theorem 3.5 (Analytic continuation of the completed zeta func- 
tion). Define the completed zeta function by 


€(s) = 59(s— 10 (3) x 8¢(s). 


Then &(s), originally defined for R(s) > 1, has an analytic continua- 
tion to an entire function and satisfies the functional equation 


é(s) =€(1— 8). 
So in full we have 


(3.4) 82D (=) C(s) = n@-D2p ($) C2), 


Proof. For the functional equation we follow Riemann’s original ar- 
gument as described in [MTB] 3.1.19]. For R(s) > 0, by definition of 
the gamma function as a definite integral, and the change of variables 


2 
co s 
Bayh 2 (4 
i gee dee 
0 


Summing over n € N, with #(s) > 1 to guarantee convergence, we 
obtain 


t = n°rx, we have 


ee) 1 
=, x? tw(a)dzx +f a? tw(a)dzx 
1 0 


or as 1 
a x? *w(ax)dzx +/ g 2 tw (=) dx 
1 1 v 


where w(z) = 7°, e-"’™@_ Note that we divided the integral into 
two pieces and changed variables by 2 +> x~/ in the second integral; 
this leads to rapidly converging integrals. The absolute convergence 
of the sum-integral justifies the rearrangement of terms. We will show 
in a moment that the function w satisfies the functional equation 


1 1 
Ww (=) — -5 — il + 21/24(z) 
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for > 0. Simple algebra then shows that for R(s) > 1 we have 


oo 
(3.5) 135 (=) (es — = +/ (228-} 4 a 39-4 u(x) de. 
2 s(s—1) 1 

Because of the rapid decay of w(x), the integral on the right con- 
verges absolutely for any s and represents an entire function of s. 
The remaining assertions of the theorem are easy consequences of the 
location of poles of 1/s(s — 1) and the invariance of the right-hand 
side of (8.5) under s ++ 1—s. It remains to verify the functional 
equation of w. For this we write 


O(x) —1 
w(x) = ( , 
with 
+00 ‘ 
O(a) = ae 


Now this series is rapidly converging for x > 0. The desired functional 
equation for w easily follows from 


O(a!) = x7 0(x), x > 0, 


which we now prove. 


Hold x constant and introduce the function f as follows: 


2,,-1 


(oo 
This is a standard Gaussian curve of the form t > e~*”. The Fourier 


transform of e~” is y > \/z/a-e-" ¥/*. With a = 1/z the Fourier 
transform of f is given by 


and, in particular, we have 


iC a eae mle 
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Now apply the Poisson summation formula [MTB| §11.4.2]. We get 
(a) =" Fn) 
=i) 


226(a) 


as required. 


The functional equation can also be expressed asymmetrically 
IMV\ p. 329]. 


Corollary 3.6. For alls £1, 
C(s) = C(1 — s)2*a*-T(1 — s) sin(s/2). 


Proof. We need the reflection principle (C.6) and the duplication 
formula (C.9) in [MV]: 


T(s)P(1—s) = 


TT 


sin 7s 
and 
I(s)l'(s + 1/2) = 72" *4T (2s). 


Applying these two formulas to the symmetrical form of the func- 
tional equation (8.4), we obtain 


¢(s) = ec 7 s\n 1/2 
1l—s 8\ sin(ms/2) s—1/2 
=r{ ; )-r(1 5) eS aaa 
~ Va2°T(1 — 8). sin(as/2) (1 — s)n°-¥/?2 


= ¢(1— s)2°x*/P(1 — s) sin(7s/2). 


We now look for zeros of the Riemann zeta function ¢(s) in the 
open half-plane #(s) < 0. If R(s) < 0, then R(1 — s) > 1. In this 
domain, the right-hand side of the functional equation is holomorphic 
and zero-free. Therefore, in the domain R(s) < 0, the left-hand side 
of the functional equation is holomorphic and zero-free. 
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The simple poles of T'(s/2) at the points s = —2,—4, —6, —8,... 
must therefore be canceled by simple zeros of the zeta function at 
these points. These are the trivial zeros of the zeta function. 


Theorem 3.7. The Riemann zeta function ¢(s) has simple zeros at 
the negative even integers given by 


Cy ee ee; ee 


These are the trivial zeros of the zeta function. 


So we know that all the zeros of the Riemann zeta function are 
trivial in the open half-plane #(s) < 0, and we have proved in Chapter 
2 that the zeta function was zero-free in the closed half-plane R(s) > 1. 

What about R(s) = 0? To find out about possible zeros on this 
vertical line, consider the functional equation with s = it,t #0. We 
obtain 


1-i 
n—*/22(it/2)C(it) = e@-VPP (=) ¢(1 — it). 
Since ¢(1 — it) 4 0, it follows that ¢(it) # 0. Therefore, the zeta 
function is zero-free on R(s) = 0. 


So we are led to ask about the zeta zeros in the critical strip 
0<R(s) <1. 

In the critical strip, the zeta zeros have the following two sym- 
metries. By the functional equation, we have 


¢(s)=0= ¢(1-s) =0. 
We already know that 
¢(s)=0 => ¢(s) =0 


where § denotes the complex conjugate of s € C. So the zeros of 
the Riemann zeta function in the critical strip come in fours, that is, 
s,8,1—s,1—%, which lie on the vertices of a symmetrically placed 
rectangle in the complex plane: symmetrically placed about the crit- 
ical line (s) = 1/2. So the Riemann hypothesis says that the width 
of this rectangle is zero; in other words, the zeros all lie on the critical 
line. 
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3.3. Some zeta values 
Lemma 3.8. ¢(0) = —1/2. 


Proof. We will use the recursive formula ['(s + 1) = sI'(s) and the 
fact that ¢(s) has a simple pole at s = 1 with residue 1. We start 
with the asymmetric form of the functional equation (8.6): 


¢(s) = C(1 — s)2°x**T(1 — s) sin(as/2). 
Multiply by 1 — s to obtain 
(1 — s)¢(s) = €(1 — s)2°x**(1 — s)I'(1 — s) sin(as/2) 
= ¢(1 — s)2°x*'T'(2 — s) sin(rs/2). 
Now let s > 1. We obtain 
—-1=¢(0)-2-T(1)-sin/2 
= 2¢(0) 


as required. 


Lemma 3.9. ¢(—1) = —4. 


Proof. To determine the value of ¢(—1), we use the functional equa- 
tion. Setting s = —1, we have 


€(—1) = €(2) 
which gives 
Va (—1/2)¢(-1) = wT (1) ¢ (2). 
We also have 
1 
1/2 
= —2T (1/2) 


reas T (1/2) 


so that 
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Lemma 3.10. ¢(—3) = 7. 


Proof. Substituting s = —3 in the functional equation 
¢(s) = C(1 — s)2°x*'T(1 — s) sin(rs/2), 
we get 


¢(—3) = ¢(4)27- 3.2 4T(4) sin(—32/2) 


ta 

— ——s+—s Qe] 
90 373 

“ke 8 

~ 90 4 

ace: 

~ 120 


Lemma 3.11. ¢(1/2) = —1.46035. 


Proof. Here, the functional equation is of no use since ¢(1/2) occurs 
on both sides. However, we can proceed as follows. Define 


n(s) =1—-2°-°43°%-4°%4.--- 
nm. ; 
n=1 


This series is convergent in 3(s) > 0. Then we have 


Gyatg aan Sf eas he a 
WL tog he fa 3s a5 | 


= ls) — cls) 
= (1-2'-*Y¢(s). 
Setting s = 1/2, we have 


2) = 5 a 


a slowly convergent series. 


Suppose that you were very naive and (like Ramanujan) knew no 
complex analysis and (like Ramanujan) had a unique and inexplicable 
mathematical intuition. You might stray from the path of rigor as 
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follows. Moving outside the domain of definition of the zeta function, 
let s = —1 in the defining equation: 


co 


Doge 


n=1 
ee) 

= 7 
n=1 


leading to the spectacular equation 


¢(—1) 


1 
ae 
This equation appeared in the first letter (received about January 31, 
1913) from Ramanujan to Hardy. 


1424+344+4--= 


This was Ramanujan’s brave, and correct, evaluation of the zeta 
function at s = —1. 


Here is another spectacular gem from the pen of Ramanujan (it 
appears in his first letter to Hardy p. 167]): 


19 + 29 +39 4.49 4.-.= 1/120. 


How would you interpret this equation in terms of the zeta func- 
tion? 
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We first need to reformulate the functional equation. Start with 
Corollary (8.6): for all s 4 1, 
¢(s) = ¢(1 — s)2°x* 1D (1 — s) sin(as/2). 
Then 
1 
1—s) = (2n)-* 
er) 
= (2n)-* Ms) sin(7s) 
sin(7s/2) 
= 2(27)~* cos(rs/2)T'(s)¢(s) 
by the reflection principle C.6}: 


T(s)P(1—s) = 


¢(s) 


sin7s— 
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The functional equation of the zeta function, established by Rie- 
mann, who was the first to treat the zeta function as a function of 
the complex variable s, is given by 


(3.6) ¢(1 — s) = 2(2m)~* cos(ms/2)T'(s)¢(s). 


Euler was led to this equation 110 years before Riemann. Now 
Euler worked with real s. For s a positive integer > 2, 1—s is a nega- 
tive integer and so the series for the corresponding zeta value diverges, 
showing that its value has to be interpreted by a summation proce- 
dure; Euler used the generating function method (Abel summation). 
Euler worked with 


n(s) =1—2°°+37%—---=(1—2'-)¢(s) 


rather than ¢(s) for better convergence and more accurate numerical 
evaluation. 


Euler, in his paper E352|1749], conjectured an equation 
which is equivalent to the functional equation of the zeta function. 
This work of Euler was forgotten, and 110 years later Riemann would 
rediscover this functional equation and provide a rigorous proof. All 
modern proofs of the functional equation involve mathematical tools 
which were unavailable to Euler, and it is remarkable that he was 
nevertheless able to predict its structure. Osler shows in his paper 
how Euler used divergent series to discover the functional equation 


for the zeta function [Osler]. 


In 1749 Euler gave a paper to the Berlin Academy entitled Remar- 
ques sur un beau rapport entre les séries des puissances tant directs 


que réciproques, E352|1749]. 


Instead of using ¢(s), Euler uses the alternating series 
co (=1)""* 
n(s) = D a 


n=1 


which is defined in the wider region #R(s) > 0. The function 7(s) is 
one simple step removed from ¢(s) as shown by the relation 


n(s) = (1 — 2°*)¢(s). 
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The functional equation (3.6) now becomes 
(3.7) (2°-1 — 1)n(1 — s) = —(2° — 1)2~* cos(ms/2)T(s)n(s). 
Euler notes the following relations E352|1749, §9]: 
1—-2+3-4+5-6+-:- _ ,1-(2-1) 
1 g2 T 32 4 vT 52 a et (2 > 1)? 
1—2?43?-4?+5%-6? + 
1 I Li = 0, 
93 T 33 43 T 53 63 
1— 23 +33 — 43453 63+... _ 1-2-3(24-1) 
1 gr 1 3m - r 54 = mrt (28 = 1)r4 
1-244 34-44+454-64+ 
1 1 1 =0, 
95 T 35 45 55 6° 
1—2°4+3°-—4°4+5°-6?4+---  -1-2-3-4-5(2®-1) 
1 96 T 36 a 56 a oe (25 = a , 
1 — 26 + 36 — 46 + 5° — 68 + 
rT I I =0, 
a7 37 Av 57 67 
1—27 497 —47 457 — 67 +... _ 1-2-3-4-5-6-7(28—1) 
1— 98 + 38 ro v 5S 3 Cie (2) =n? 
1— 2% +38 —48 458-68 + 
1 I ais = 0, 
29 T 39 a9 ' 59 69 | 


1— 29 + 39 — 49 +. 59 — 69+... 
1 sit gio Th + sto gio men 
-2-3-4-5-6-7-8-9(2'° — 1) 
(29 — 1)! 


g(l—n) _ (1) (2" — 1) — I)! if nis even 
a) Oia 


=0 if nis odd. 


These equations are now rewritten in the form 


nl—n) -(a— 12" =1) ™m 
may (QP Tee 
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Euler, in his article E352|1749, §10], now writes 


‘Je hazarderai la conjecture suivante, quelquefoit l’exposant n”: 


ml—n)  —l+2+3-9(— DQ" —1) ™ 
nn) ("1 = I)a = 


With s = n, this is precisely (3.7)! 


Euler verified his conjecture for n = 1/2, which is 


1- Lo + a _ 1/41/2 fee —T(1/2)(21/2 =i) o 
1—1/21/2 41/31/2 1/41/24... (Q-1/2 — 1) 71/2 cos 7 


This is the only case considered by Euler in which both the series in 
the numerator and the denominator converge. He used ['(1/2) = ./7, 
and the right-hand side immediately reduces to 1. Euler remarked 
that his success in this case makes his conjecture very convincing. 


See E352|1749, §14]: 

“There should not be any further doubt about our conjecture, 
having verified it not only for all the cases where the variable n is 
an integer, be it n is positive, or n is negative, but also for the case 
n= 1/2.” 

He goes on to test his conjecture for n = (2k + 1)/2 with k EN. 
For n = 3/2 he obtains 


1-2 4 + bre 34+ V2 
- Bae le oes eet _ = ms = 0.4967738: 
a/i B/8 4a Bs eve | ORV? 
see E352|1749, §15]. With n = 15/2 he gets 
1 — 2/34 99V3— 48/44... 1-3-5-7-9-11-13(128V2 — 1). 
1 1 1 — ’ 
tare sya arya 27(128— v/2)n" 


see E352|1749, §16]. 
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3.5. The Euler constant revisited 


Euler discovered the constant y, defined by 


This is a fundamental constant, now called Euler’s constant. We give 
it to 50 decimal places as 


7 =0.57721 56649 01532 86060 65120 90082 40243 10421 59335 93992 - - - 


There is a whole book devoted to the Euler constant, see Julian 
Havil [H). 


The Euler constant y is equal to the constant term of ¢(s) at 
s = 1, as we now proceed to show. 


Theorem 3.12. In the Laurent expansion of ¢(s) about s = 1, we 
have 


((s) = 7 +74 O(8-1). 


Proof. First we note that 
fs -[{% [7 2Qdx | +f (n — 2)dx f (n —1)dx 
a" : n—-2 a? : n-1 a? 
Pal Ale he 
= +.+-4 
Udy tlo v In-2 v 
1 1 1 1 
=1 + 2 4 
3+2(5-3) +8 (5 il Gog) 
v(a-3) 
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Now we apply Theorem [I.10] We obtain 


wy (0o-ph)=1- fe 


In the next result, we expand the gamma function I'(s) as a Taylor 
series about s = 1 and show that the coefficient of the linear term is 


~+4. 
Theorem 3.13. We have 


I'(s) =1-—7(s—1)+ O(s—1)?. 


Proof. This is the Taylor series expansion for '(s) around s = 1: 
I(s) =T(1) +I’ (1)(s — 1) + O(s — 1)?. 
We have ['(1) = 1. In addition, we need Euler’s formula 
OO aaa 


see E368|1765 §14]. To prove Euler’s formula, we can proceed 
as follows. Introduce the digamma function 


W(s) := a = £ logs). 
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The digamma function is meromorphic and satisfies the difference 
equation 


U(s 1) =4(s) += 


which is inherited from the recurrence relation ['(s + 1) = sI'(s). A 
simple iteration now yields the equation 


n-1 
1 
va) = v(e+n)- 
j=0 
from which follows 
n-1 1 
v(2)= im, (nen 


This limit is valid for all ¢ € C\Z<o, using the fact that w(x) —log x 
goes unformly to 0 in the right half-plane as §t(a) > oo; see [L} 3.1]. 
This formula gives w(1) = —7y as required. 


The Euler constant is inextricably entwined with the zeta function 
and the gamma function. Here’s another example, which relates 7 to 
the zeta values at the natural numbers greater than or equal to 2. 
The result is due to Euler E583/1776, 86]. 


Theorem 3.14. 


ery 
n=2 


Proof. Recall the definition of the Euler constant y which was given 
by 


a 
7= Jim & een). 


m=1 
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We therefore have 


k=2 

[oe) foe) 1 
aia ep Der 

k=2m=2 

ie a ol 
ae De a 

k=2 m=2 

at (hy = 
=1-\° ) 

k=2 


which is the required result. 


Select a prime number p. Given x € Q, we have its unique fac- 
torization as powers of prime numbers, « = p* - pe ---pkn where 
DP, P1,-++;Pn are distinct primes and k, k1,ko,...,k, € Z. The p-adic 
norm is defined as follows: 


c=p*-ph...pkm => |2|,=p 


—k 
The p-adic norm makes a beeline for the exponent of p and then 
reverses it. The completion of Q with respect to the norm | - |, is Q,, 


the field of p-adic numbers. 

We also have the absolute value |a| which will presently be de- 
noted |x|... The completion of Q with respect to the absolute value 
|-|oo is the real field R. 

We recall that a place of Q is a completion of Q. We have the 
finite places Q, and the infinite place R. 

Euler’s constant y appears in the Laurent expansion of the Rie- 
mann zeta function at s = 1; here the zeta function is a function 
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encoding the properties of the finite primes. Euler’s constant also ap- 
pears in the Taylor expansion of the gamma function at s = 1; here 
the gamma function encodes properties of the infinite place. 


The Euler constant occurs in two apparently different contexts: 
one associated to the finite places of Q and the other associated to the 
infinite place of Q. There is a well-known connection between these 
two contexts, given by the product formula for Q over all places, finite 


and infinite: 
II zip = 1. 
2<p<oo 
The functional equation can be written 


C*(s) =¢*(1—s) 
where 


C"(s): =P (Z) c(s) 


_ sp (8 1 ee 
=r (I(t) 
If we write 
G=(i-7 7) 
carver(), 


ll] -e¢-. 


2<p<oo 


This point of view is the point of view of Tate’s thesis [Tate], 
and indeed of Weil’s book [We?2]}. 


then we have 


3.5.1. The Euler constant and periods. This section is inspired 
by the article of Kontsevich and Zagier and by §3.14 in the 
article by Lagarias [L]. Numbers which are not algebraic are called 
transcendental. There is a huge difference in size between alge- 
braic and transcendental numbers: the set Q of algebraic numbers is 
countable and the set of transcendental numbers is uncountable. 


There is, however, one further important class of numbers, lying 
between Q and C, which is missing in the above classification. This 
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“new” class of numbers, the periods, seems to be the next most 
important class in the hierarchy of numbers according to their arith- 
metic properties. The periods form a countable class. Periods appear 
surprisingly often in various formulas and conjectures in mathematics. 


Definition 3.15 ({KZ} §1.1]). A period is a complex number whose 
real and imaginary parts are values of absolutely convergent integrals 
of rational functions with coefficients in Q, over domains in R” given 
by polynomial inequalities with coefficients in Q. 


We will denote the set of periods by P. It is obviously countable. 
Sums and products of periods remain periods, so the periods form 
a ring. Maxim Kontsevich and Don Zagier gave a survey of 
periods and introduced some conjectures about them. 


The irrational number V2 is a period, for we have 


V2= dx. 


Qa2<1 


The periods are intended to bridge the gap between the algebraic 
numbers and the transcendental numbers. The class of algebraic num- 
bers is too narrow to include many common mathematical constants, 
while the set of transcendental numbers is not countable. Periods are 
usually transcendental numbers, but they are described by, and con- 
tain, only a finite amount of information and (at least conjecturally) 
can be identified in an algorithmic way. 

The set of all periods forms a ring P which includes the field Q 
of all algebraic numbers. It follows that if Euler’s constant is not a 
period, then it must be a transcendental number. We can express 7 
as a period by any of the following integrals: 


: 7 dz ~~ dx 
_ _ /| _ 72d — = 
T= // dady = 2 [ 1— 2x7dx Loe [= 


u2+y2<1 


or also, after multiplication by the algebraic number 22, by the contour 


integral 
d 
271 = ¢ os 
z 


in the complex plane around the point 0. 
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Two other famous numbers which have special notations are 

1 n 

e= lim (1 + =| = 2.7182818..., 

noo n 

the basis of natural logarithms, and Euler’s constant 

1 1 
y= lim (14 5 b+ 2 togn) = 0.5772156... 

nN—-0o 2 n 


but these two numbers (conjecturally) are not periods. It is known 
only that e is transcendental. 


However, there are many examples of periods besides 7 and the 
algebraic numbers. For example, logarithms of algebraic numbers are 


periods; e.g., 
2 dx 
log2= | —. 
qe 


Many infinite sums of elementary expressions are periods. For exam- 


ple 
1 1 
¢(3) = 1+ 55 + gg +--+ = 1.2020569... 
has the representation 
dadydz 
3 — —— 
(3) Teor 


0<r<y<z<l 


as an integral, and more generally, all values of the Riemann zeta 


(3) = 


n>1 


function 


at integers s > 2 are periods, as are the “multiple zeta values” 


1 
G (Sigetay Sh) = net... nse 
O<ny<e-<ng k 
with s; EN, s, > 2. 


Values of the gamma function 


Co 
I(s) = i i? eae 
0 
at rational values of the argument s are closely related to periods: 


I(p/q)t <P 
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with p,q € N. For instance, 
T(1/2)? =" 


and 


7 dx 
T(1/3 a atlSgl/ta | = 
oa riers 


We have, in this book, met ¢(2) on several occasions. Here is a 
realization of ¢(2) as a period over Q; see [HMS| Proposition 12.1.7]: 


, J 60}, 


epeyer (L-—2)y 


In 2001, Kontsevich and Zagier posed the following problem : 
exhibit at least one number which does not belong to P. 


Of course such numbers exist since P is countable. Solving this 
open problem would be the analogue of Liouville’s achievement in 
the 19th century when he constructed the first explicit example of 
a number which could be proved to be transcendental. Even more 
desirable would be to emulate the achievements of Hermite and Lin- 
demann and prove that some specific numbers of interest, like e or 
1/x, do not belong to P. 


In 2008, M. Yoshinaga wrote down a nonperiod 
a = 0.77766444.... 


This is a computable real number which is not a period. For a recent 
account of this result, see [HMS| §16.6]. 

There are many famous unsolved problems about the nature of 
the Euler constant. The most well-known one is 


Conjecture 3.16. Euler’s constant is irrational. 


This is a longstanding open problem. A recent and stronger ver- 
sion of it is the following. 


Conjecture 3.17. Euler’s constant is not a period. In particular, 
Euler’s constant is transcendental. 
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Definition 3.18. An exponential period is an absolutely convergent 
integral of the product of an algebraic function with the exponent of 
an algebraic function, over a real semialgebraic set, where all polyno- 
mials entering the definition have algebraic coefficients. 


Here is the classic example of an exponential period: 
Jr = exp(—2”)dz. 
Recall the definition (8.2) of the gamma function for R(s) > 1 
I(s) = [ er de 
0 


With a real variable x this becomes 


Ii(z)= , as aad 
0 
which is precisely Euler’s formula E368|1769 $11]. 


Differentiation under the integral sign yields 
(j= | e *logadz. 
0 
Combining this with Euler’s result E368|1769 §14] 
(1) =~; 


=| e “ logadz. 
0 


Substituting —log a = f. y ‘dy for0<a<1and 


log = f y 'dy 
1 


for x > 1 yields the formula (due to Kontsevich): 


ve ff Sane [0 [ae 


This elegant formula ne that 7 is an exponential period. 


we have 


Chapter 4 


The explicit formulas in 
analytic number theory 


In the present context, an explicit formula is an equation of the fol- 
lowing form: 


a counting function depending only on the prime numbers 


=a formula depending only on the zeta zeros. 


Such an explicit formula is said to exhibit the reciprocity between 
the prime numbers and the zeta zeros. 


The first such explicit formula was discovered by Riemann in his 
memoir of 1859 [R]. Later, von Mangoldt proved his own version of 
the Riemann formula. It is customary to focus attention on the von 
Mangoldt formula, but in this book the focus will be on the Riemann 
formula, which deals explicitly with the difference r(x) — Li(x). 


On the right-hand side of an explicit formula we have the set 
of zeta zeros. For reasons which will become clear in this chapter, 
the imaginary parts of the nontrivial zeta zeros are of particular im- 
portance. We recall that the zeta zeros occur in complex conjugate 
pairs. For those zeta zeros p with S(p) > 0, it is customary to write 
the imaginary parts as 71 < yo < 73 <-°:- 

In the book by Mazur and Stein [MS], the set {71, 72, 73,..-} is 
called the Riemann spectrum of the primes. 


57 
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In this chapter, we will state and prove the two explicit formulas 
and discuss their relationship with the Riemann hypothesis. The 
Riemann spectrum may be interpreted as a sound spectrum, and we 
investigate how this sound varies with time. 


4.1. The von Mangoldt explicit formula 


We recall Definition [2.5] of the function A (n) where n € Z. We have 
A(n) = logp if n = p* and A(n) = 0 otherwise, where p is a prime 
number. Recall also the logarithmic derivative of ¢(s): 


Theorem 4.1. The von Mangoldt explicit formula is given by 


(4.1) v(a)=2 >= 5 log (1 =) log 2a 


p 


where p is the zeta zeros which lie in the critical strip 0 < R(s) < 1. 


Note that the left-hand side of the von Mangoldt explicit formula 
depends only upon the prime numbers, and the right-hand side 
depends only on the nontrivial zeros p of the Riemann zeta function. 
We can, if we so choose, incorporate the trivial zeros into this formula 
in the following way: 


—n 


1 ae, 
«(1-3)-E5 


n=1 
1 1 ae ie 
as | 1—-—)= 
—- 2 os ( =) d 2n 
co —2n 
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The integers {—2n : n = 1,2,3,...} are the trivial zeros of the Rie- 
mann zeta function. So we can now write the von Mangoldt explicit 
formula in the following way: 


or 


where the summation is now over all the zeta zeros. 


At this point, we are strongly influenced by $10.2 in Jeffrey Stop- 
ple’s book [St]. This reference contains some excellent diagrams, 
especially Figure 10.2. 


The issue for us at the moment is to seek a simple geometrical 
interpretation of the von Mangoldt formula. Recall that the nontrivial 
zeta zeros occur in complex conjugate pairs p and p. If p = 6+ 17, 
then noting that 


ot = ety loge and |o\? = oi fen 


we have 


aP aP ghtiy gh ty 


pp Pry pay 
po. B= ine + (B+ ie 
B? +9 
af i -i (apd -i 
= gaa (Ble ta '?) —ig(a'? — a 7)) 
2xP ; 
= Brae (6 cos(y log x) + ysin(y log x)) 
2x8 
=a G cos(7 log x) + ue sin(y log »)) 
lel \ Tel |p| 
QF 
- a cos(y log x + @) 


where @ is such that cos@ = '/ |p| and —sin@ = ¥/ |p]. 
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We could use a logarithmic scale. Let t = log; that is, « = e’. 
Then we have 
ere = FE cea Gy: 
p Pp |p 
This can be interpreted as an oscillatory term, with t as time variable: 
~ the amplitude is 2e°'/|p|; 

— the phase shift is 6; 

— the period is 27/7; 

— the frequency is 7. 

Each complez conjugate pair of nontrivial zeta zeros contributes 
an oscillatory term, with a definite amplitude, phase shift, period, and 
frequency. 

For those zeta zeros p in the critical strip with S(p) > 0, we will 
denote the imaginary parts by y1 < y2 < 73 <-:+: <n <--+:. So we 
have pp = By +i yn. It is then clear that as n — oo we have 


~the nth amplitude 2e%""/|p,,| > 0; 
— the nth phase shift 0, = cos~!(8n/|Pn|) > 7/2; 
— the nth frequency yn, — oo. 


We list the first few nontrivial zeros of the zeta function: 


1 © 1/2 + (14.12473) i, 
pz © 1/2 + (21.02204) i, 


3 © 1/2 + (30.42488) i, 


pioo 1/2 + (236.52423) i, 


The case when 8 = 1/2. Suppose for the moment that the Rie- 
mann hypothesis is true. Then all the nontrivial zeta zeros lie on the 
line R(s) = 1/2. In other words, we have that G= 1/2. So we can 
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now write / 
Pp Pp Qei/2 
aula ea eens cos(y log x + 6). 
pp lp| 
The envelope of this curve is then the parabola y = +2x1/?/ |p| or 
2 
y? = 4a/ |p|". 


As |p| increases, we can see that the parabola y? = 42/ |p|” gets 
closer and closer to the z-axis. 


For a vivid graphical display of the von Mangoldt formula, go 
to www.demonstrations.wolfram.com and enter Using Zeta Zeros 
to Compute the Chebyshev Psi Function in the search tab. If you 
download the Wolfram Player app, then you can vary the parameters. 
Even with the first 200 pairs of zeta zeros, you can see how well the 
von Mangoldt formula hugs the step function w(z). 


4.2. Can you hear the Riemann hypothesis? 


We consider once again the von Mangoldt explicit formula 


p 
v(e)=2- D> +40). 
¢(p)=0 e 

The left-hand side of this equation is a step function. On the right- 
hand side, we have already shown that each complex conjugate pair of 
zeta zeros gives rise to an oscillatory term, with amplitude 2e°'/ |p]. 
If we fix a zeta zero p, then this amplitude is increasing as ¢ increases. 
If we think of ourselves as listening to a sound wave of this amplitude, 
then we will experience a crescendo; i.e., the sound will get louder. 
We attempt to quantify this phenomenon. An a-crescendo is the ratio 
a of the amplitudes taken over a time interval of one second. With 
initial time T, the amplitude ratio is 


=) / (=) (=) ( |p| ) 8 
— =e 
lp| lp| lpl 2ehT 


which is independent of the initial time T. Each conjugate pair of 


nontrivial zeta zeros gives rise to a pure tone with an e?-crescendo. 
So if the Riemann hypothesis is true, then each conjugate pair of 
zeta zeros gives rise to a pure note with a ,/e-crescendo, and as t 
increases the loudness of this note increases at a constant rate. When 
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we take the superposition of all the zeta zeros, we obtain a chord. We 
introduce the term chord of the zeta zeros. The first nontrivial zeta 
zero gives us the fundamental frequency of this chord; the harmonics 
are given by the other zeros. 


Consider NH (Negation of the Riemann hypothesis). In particular 
suppose that there exists a single zeta zero sq which is off the critical 
line: 

80 = Bo + 20, Bo #1/2. 
Let 89 = s +6 with 5 4 0. For this pure note, we have a \/é - e°- 
crescendo. 


If 6 > 0, then this pure tone will get louder more rapidly than all 
the other pure tones in the chord of the zeta zeros. A listener with 
very acute hearing should be able to detect this! 


If, on the other hand, 6 < 0, then this pure tone will get louder 
more slowly than all the other pure tones in the chord of the zeta 
zeros. All the other tones are performing a crescendo at exactly the 
same rate /e. 


On the relation between loudness and amplitude. We often express 
amplitudes in logarithmic units called decibels. If a is the amplitude 
of a signal, then we can define the decibel (dB) level d as 


d = 20- logy) a/ao 


20 - logyg a — 20 - log;, ao 
where ao is a reference amplitude. A convenient choice of ag is 
ap = 107°. 
When the amplitude a is multiplied, as in the present context, by e®, 
the decibel level rises by 
20 - logy, e® = 20 - Blog, e © 206(0.4343) = 8.6868. 


When 6 = 1/2, the rise in decibel level is ~ 4.343. Amplitude 
is related in an inexact way to the perceived loudness of sound. But 
amplifying a signal by 3 dB, say, will fairly reliably make it sound 
about one “step” louder. So, in our case, amplifying the signal by 
4.343 will certainly be audible. For more information on the science 
of sound, in particular the perception of sound, see [RM WY} p. 99]. 
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If the Riemann hypothesis is true, then we will obtain a steady 
crescendo for each pure tone in the chord of the zeta zeros: for each 
unit of time, the dB level will rise by by 4.343. 


4.3. Comparison with Fourier series 


We turn once again to the von Mangoldt formula: 


v@)=2-) = +00) 


p 


where the summation is over all the zeta zeros. If the Riemann hy- 
pothesis is true, then we can write 


v(e)=a-)) 


ll2+iy 


Va rue 
We know from above that 
xP P 2x? 
+—= B cos (ylogx) + ysin (ylog2)). 
rt = Gey (B08 (ylog.2) +7 sin (ylog2)) 
Now 
So sin ylog x (+- ae giz)- S- sin ( (ylog x) ee Cai ie ails 
or) B a ne; epee ty) 
2 
= So sin( (ylog x) \(- ae ) 
a; (B29?) 
= O(1) 
since 6? <1 and eee 1/7? < oo. Similarly, we have 
] 
ye ee) og x) ~O(1). 
aa? err 


We infer that 


ae) = ee sin (¥ log x) 
ie = 22, 3 + O(1). 
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If, as before, we write the zeta zeros p in the critical strip with S(p) > 
Oas $+in,s +in2,5 + 93,-.. with 0 < y < 72 <43 <--:, then 
we have 


v(a)-—2 sin (Yn log x) 
Vr d, Yn ) 
If we set t = logx, then this may be written 
y(a)—« <= sin(ynt) 
4.2 =-2 + O(1). 
(4.2) me d - (1) 


Now we have the following standard Fourier series: 


2a—1 \ sin (27a) 
=—2 
2 d 27 

with x € [0,1]. We may extend the domain of x to the whole real line 
as follows: 


2a — 2 [a] sin ( (20na) 
4. =-2 
oo 2 > 27m 


At this point, it is interesting to compare with (42). In the 
case of we have the frequencies 27, 47,67,.... These frequen- 
cies are all integer multiples of a single frequency (the fundamental 
frequency). In the case of we have the frequencies 71, Y2, Y3,---- 
Here the situation is quite different. It is conjectured that these fre- 
quencies are all linearly independent. 


Even so, the superficial similarity of (43) and (4.2) has led us to 
interpret the right-hand side of (42) as a sound wave, a superposition 
of pure tones. Such a sound file can be made with the Mathematica 
package called Audio. 


Historical comment. The equation was the starting point in 
Littlewood’s proof of his famous theorem concerning the sign oscil- 
lation of 7(a) — Li(#). This equation appears as equation (5) in his 
very terse Comptes Rendus Note [Li], published shortly before the 
outbreak of WW1. He writes 7 = log x so his equation (5) appears as 


ae v_ 2 ety O(1). 
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4.4. Proof of the von Mangoldt formula 


Where does the von Mangoldt formula come from? Recall the defini- 
tion of the Euler product, which was defined in Chapter 2 as 


Il (1 = a ' 


Pp 


This is a product over primes. There is a totally different way of 
thinking of the zeta function, which is to think of it as a polynomial 
of infinite degree. In that case, leaving aside its simple pole at s = 1, 
it should be possible to split it into its linear factors. This heuristic 
argument at least makes it plausible that, for R(s) > 0, the correct 
infinite product (the Hadamard product) is 


seca (2) Tse 


where the product is taken over all the zeta zeros p. You can see the 
simple pole at s = 1 and each zero at s = p. 


We have already met the Weierstrass factorization theorem in 
Section 3.1. It’s time to state it explicitly p. 220}. 


Theorem 4.2. Let f be an entire function with {21, 22,...,2;,.-.} 
as its zero set, each z; being counted as often as its multiplicity. Let 
m be the multiplicity of 0 (m may be 0). Then there exist integers 


P1;P2,+-+,pj,--. and an entire function g such that 
f(z) = exp(g ™.J[E pj (2/25) 
j=l 
where 
Eo(z) =1- 
E\(z) = aoe 


E,(z) = (1—z)exp(z + 27/2 +--+ +2?/p). 
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The Hadamard product fits into this theorem by taking 


f(s) = (6 — 1)*4(s), R(s) > 0, 
m= 0, 

g(s) = slog(27/e) — log 2, 
pj = 1, 


for all 7 = 1,2,3,.... 


The next crucial step is to equate these two infinite products for 


C(s): 
(444) [JQ-p*) t= au (=) [] G-s/p)e%"”. 


P 
Note that all the primes are on the left-hand side and all the zeta 
zeros are on the right-hand side. The explicit formulas in number 
theory are all obtained from (£4) by various clever manipulations. 
In the course of these manipulations, the duality between primes and 
zeta zeros is preserved. 


We will now multiply through by the pole at s = 1: 


1 (2n\* 
(4.5) (s-1)¢)=5 (=) []a-s/oer. 
p 
We now take the logarithmic derivative of (4.5); i.e., we first take the 
logarithm of each side and then differentiate. Spelling out the details, 
we get, for the left-hand side, 


G(s) + (s = 1)¢'(s)) fp 


(s—1)¢(s) ss -1 Gs) 
If f (s) = (1—s/p) e*/?, then 


i) (pel? + L—s/oyer "Ja 
f(s) (1—s/p) e8/e 
—s/p" 
~ 1-s/p 


p(s—p) 
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If f(s) = 4 (22)°, then we get 


= log(27) — 1. 


Putting all of this together, we have that the logarithmic derivative 
of the Hadamard product is 


1 0s) ys ae 
oy oa 1. 


Then we we have 


(s) __s ee 
@ i op 
By Theorem [2.6] we have 
C(s)_ yn An) 
(4.6) Oe d a 
We have 
“A(n) os 8 
S- Sor ee Die eer, + log (277) . 


Note that the fundamental duality between primes and zeta zeros 
is preserved in this equation: the left-hand side depends entirely on 
the primes, and the right-hand side depends entirely on the zeta zeros. 
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We get to work on (46): 


n=l ne n=1 n® 
es 1 1 
= Dae Gs aa) 


I 
| 
_— 
= 
— 
a 
i 
i 
a 
8 


II 

| 

wn 
_— 
So 
~ 

8 
o 
i 
i 

8 


where M denotes the Mellin integral transform, which is defined as 
follows. 


Definition 4.3. The Mellin transform (Mf) (s) is defined as 


(Mp)(8)= fo f(e)a- de 
1 
under certain conditions on the function f and the variable s. 


The Mellin transform has the following properties: M is linear 
and injective. 

We now take the Mellin transform of the right-hand side of the 
von Mangoldt explicit formula. That is, we take the Mellin transform 
of 


Pp 
2 — Y= — log (27). 


p 


We require an elementary property of the Mellin transform. 
Namely, 
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for R(s) > R(a); see Exercise 5.3 in the final section of this book. 
Substituting a = 1, p,0 in this formula, we get 


—s(M (2)) (8) = —, 


=e a 
—sM(— log (27)) (s) = log (27) . 


Putting all this together, we have that 


—sM (- > = ~tx(2)) (s) = i ° : +0 a= + log (27) 


p 


We infer that 


Now we exploit the fact that the Mellin transform is injective and 
conclude that 


xP 
z— >) — —log(2m) = o(2), 
Pa 
which is the required result. This completes the proof of the von 
Mangoldt explicit formula. 


The von Mangoldt explicit formula comes from the inverse Mellin 
transform of the logarithmic derivative of the Hadamard product for 
the zeta function. 


4.5. The logarithmic integral Li(z) 


As preparation for the Riemann explicit formula, we need a careful 
discussion of the logarithmic integral Li(z) of a complex variable. It 
is, of course, astonishing that Riemann used the function Li(z) in the 
very early days of complex analysis. Riemann does not define Li(z) 
in his article of 1859, so it was left to later researchers to figure out 
the definition. We start with the exponential integral as in Lebedev 


§3.1]. 
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Definition 4.4. Define the exponential integral as follows: 
z ew 
Ei(z) = —du, |arg(—z)| <7, 
oo W 
where the integration is along any path L in the w-plane with a cut 
along the positive real axis. 


The function Ei(z) is holomorphic in the plane with a cut along 
the positive real axis. A possible choice of the path of integration is 
the infinite line segment 

—00 < R(w) < R(z), Sw) = S(z), 
passing through the point z and parallel to the real axis. 

According to p. 32], the values of Ei(z) on the upper and 
lower edges of the cut are, respectively, 

Ei(u + i0) = Ei, (u) — ix, u > 0, 
Ei(u — 20) = Ei, (w) + iz, u > 0, 
where Ei;(u) is the real function defined by 


y” 


Tin’ u > 0. 


Eiy(u) = 7+ logut+ » 
n=1 


This is in preparation for the next definition. 


Definition 4.5. Let s = u+it. Define the logarithmic integral as 
follows: 


Li(x*) = Ei(slogr)+im if S(s)>0 


S$ 
S$ 


=Ei(slogr)—in if S(s) <0 


On the upper and lower edges of the cut we have 
Li(z"*"°) = Ei((u + 20) log x) + im 
= Ei, (wlog x), 
Li(a"—*°) = Ri( 
= Fi, 


(u — 20) log x) — ix 

(ulog x), 

which effects a smooth transition across the cut for the function s > 
Li(a*). 


4.5. The logarithmic integral Li(z) 71 


Note that Definition [4.5] is not the definition of Lebedev 
§3.4]. However, it appears to be the definition favored by Riemann: 
on this point, see p. 30 of the book by Edwards and the note 
on p. 186 of [R]. Definition [4.5] is also the definition favored by the 
Mathematica programmers. 


Let x > 2. We now consider the curve 
tr—> Li(at/2+#) 


whose domain is the critical line #(s) = 1/2. Definition [4.5] allows us 
to analyse this curve. It will turn out to be a double spiral. 


Theorem 4.6. The map t+ > Li (ae) defines a double spiral for 
which 


jim Li Coe =in and jim Li Cog =-1T. 
—>0o 00 


The upper arm spirals counterclockwise to the limit point it and the 
lower arm spirals clockwise to the limit point —ix. The upper arm is 
the mirror image, with respect to the real axis given by S(s) = 0, of 
the lower arm. The double spiral meets the real axis at the point 


(1 
Ei; ((log x) /2) = 7 — log 2 + log log x + 3 dosa ? 
2"nIn 
Proof. Fix t € R. Let s = (u+it)logaz. We have 
ds __— du 
s urtit 
and —oo < u < 1/2. We have 
5+it) log es 
(4.7) Ei((1/2 + it) log x) i —ds 


(4.8) 


3 UW 
(4.9) =_" We ~ du 
oo UTI 


ll 
— 
Q nie 
S on 
+ a, 
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and 


|Ei((u + it) log x)| = 


IA 
Je 
, 
eas 
bo 
8 
iS 
Q 
S 


1 gil? 
 t (Z) 
Then Li(a/?+) —+ in as t —> oo. Now suppose that s = 


(v — it) loga. Then again we have 


z UV m—it 
2 


Ei((1/2 — it) logx) = / dv 


ge 0 St 


= Ei((1/2 + it) log x) 


and 


|Ei((1/2 — it) log x)| = 


1 gi/2 
=F (F) 
Therefore we have Li(a!/?-“) —+ —im as t —> oo. 
The phase of the complex number Li(x!/?+“*) is, in view of (£7), 
dominated by the phase of x”, which is tlogz. This phase steadily 


increases as t —>+ oo and indicates that we do indeed have a spiral 
curve. 


The zeta zeros which lie on the critical line can now be plot- 
ted on the double spiral. For a vivid graphical display, go to www. 
demonstrations.wolfram.com and enter Logarithmic Integral on the 
critical line in the search tab. If you download the Wolfram Player 
app, then you can vary the parameter x and the range of t. You can 
also zoom towards the limit point ix. Strongly recommended! 
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Consider the complex conjugate zeta zeros p,p. Let p = 8 +77. 
We have 
B gig tt 


dv 


Ei((8 — it) logx) = i 


hk O40 


= Ei((6 + it) log x) 
and then 
Li(x?) = Li(a?) 
which is compatible with the reflection principle. We also have 
Li(x”) + Li(x?) = 2R(Li(x?)). 


It is these real numbers, one for each complex zeta zero p, which 
make their contributions to the right-hand side of the Riemann for- 
mula, as we will see in the next section. 


4.6. The Riemann formula 
Definition 4.7. Let z > 2. We define 


1 g/k 
Ore eae: 
k>0 
where x (a) denotes the number of primes less than or equal to «. 


This is a finite series: the series terminates as soon as x!/* < 2. 


Theorem 4.8. The Riemann explicit formula is given by 


J (x) = Li(x) — Soli (2°) — log 24 i t(2 S log t 


where the summation is over the complex zeta zeros p. 


Before moving on to a proof this formula, we attempt to shed 
some light on the individual terms, sometimes stating the obvious. 
e The first term Li(«) is monotonic increasing. 
e The second term is the most interesting. When we pair the 
terms coming from complex conjugate zeta zeros, we obtain the 
conditionally convergent series 


-25" R(Li(x*)) 
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summed over the zeros in the critical strip. Is the sum of this 
series positive or negative? Does it dominate, or is it dominated 
by, the first term Li(x)? We will address this issue shortly. It is 
the key to understanding this formula. 

The third term is the constant — log 2. 

The fourth term is 


as dt 
ae ie t (t2 — 1) logt’ 


This term f(x) is monotonic decreasing, with f(2) = 0.14001 
and 


Jim f(x) =0. 


The fourth term has an interesting interpretation. If we use the 


substitution u = logt, then we can write 


With a second substitution t = e 


roe dt Sa aa du 
i t (t2 — 1) logt =[ wecD 
Py his eX du 
-[ 1—e724%@ y 


fo <a 
io € 2nu dy, 
log x U 


n=1 
oe) 


3 i edu 
log x U 


n=1 


—2nu we obtain 


oof en 2mUdy Sf? dt 
a ae 
ae 


n=1 


which is the second term in the Riemann formula, but now the sum- 


mation is over the trivial zeta zeros. So we can, if we choose, write 


4.6. The Riemann formula 75 


the Riemann explicit formula in the following concise way: 
J (x) = Li(a) — $> Li (2?) — log2 
p 


where the sum is over all the zeta zeros. 


We now have to extricate 7(x) from the finite series J(x). The 
method is Mobius inversion. 


Definition 4.9. We define the Mobius function yu as follows: 


(1) wl) =1. 
(2) uw (n) = 0 if n has a squared factor. 
(3) Given distinct primes p1,...,p%, then 


p(p1 - p2-+* Pe) = (—1)*. 


Lemma 4.10. Let yw denote the Mobius function. Given a positive 
integer n, we have 


Su (d) =0 
d\n 


where d\n means that d divides n. 


Proof. Suppose that n > 1 and that n = p{'p5?---pZ*, where 
P1,-++;Pk are distinct prime numbers. Then we have that 


k 
TH =14 Dao) + EE ato) +-+4([Lo] 


d\n i>1 i>1 j>1 
k 
=1-k+4 +++++(-1)' 
(B) tt 
=(1-1 
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We immediately exploit this result. Let k be determined by the 


conditions z,2!/2,...,a!/* > 2 but a!/(*+) <2. Then we have 
k k k 1/ij 
3 (3) (21/) = s H3) + 1 (2! 3) 
j=l J g=l J w=1 ‘ 
k? 
= u(d) (/") 
n 
n=1 d\n 
= 1(2) 


with the change of variables n = ij. This expression for m(a) is a 
finite series in the terms J(a),..., J(«!/"). Combining this with the 
Riemann formula, we obtain 


t? — 1) logt 


J (2) = Li(e) ~ J Li(a?) ~ tog2 +f mt 


where the sum is over the complex zeta zeros. 


This is one of the great formulas of all time, expressing as it does 
the number of primes up to x in terms of the zeta zeros. 


We recall the classical theorem: 


This result is stated by Euler in his book Chapter XV, p. 235] 
and was subsequently proved by von Mangoldt; on this point, see the 
discussion by Edwards p. 92]. 

Thanks to this result, the constant —log2 will disappear after 
Mobius inversion. It is customary to ignore the fourth term in the 
Riemann formula, the very small term which tends to 0. In that case, 
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we have 
k 
(4.10) n(x) = S> pin) Li(x'/") — 25° R(Li(2?/)) 


n=1 p 
where the summation is over the zeta zeros in the upper half-plane. 


We shall look carefully at this extraordinary formula for large 
values of x. It certainly looks as if Li(a) is going to win; i.e., it looks 
as if r(a) < Li(a) for all x. For Li(z) is immediately followed by three 
large negative terms corresponding to the primes n = 2,3, 5. 


The only way that 7(«) is going to beat Li(z) is if the zeta zeros 
arrange themselves on finitely many double spirals S,,...,S, in a very 
particular way. According to Littlewood’s theorem, this is precisely 
what can happen! From a naive point of view, one would expect a 
great deal of cancellation, and indeed most of the time, this kind of 
cancellation does take place. But there are certain intervals on the 
real line for which the real values on these spirals reinforce each other 
sufficiently to dominate Li(x). This reinforcement must be negative 
in order to dominate the positive term Li(zx). 


Such intervals occur only when x is very large, as we will see in 
the next sections. 


Choose and fix x > 2 and then think of x as large. To be clear 
about the finitely many double spirals, first note that «°/I = (a!/1)¢. 
Now consider the double spiral 


S:tt> Li (y3+*) 


and allow y to take the values az, #!/2,21/3,...,2!/*. This will create 


double spirals S1,...,Sx. 

These double spirals are getting smaller, for the real numbers 

Sy, M real axis 
decrease as n increases; see Theorem (46). To be precise, we have 
Sp, Areal axis = Ei, ((log x!/”)/2) = Ei, ((log x) /2n). 

B(9) 
numbers Li(a?/ ”) will lie on the double spiral S,, ai emies pisa zeta 
zero on the critical line. 


We also have to scale each double spiral S; by . The complex 
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Here is a thought experiment. Imagine looking at a large number 
of these double spirals S;,...,S,. Assume the Riemann hypothesis, 
and plot the complex zeta zeros on all of these spirals. If you look at 
this configuration from a great distance, then these millions of points 
will usually arrange themselves equally on each side of the imaginary 
axis S(s) = 0. However, for certain values of x, these millions of 
points will have a tendency to arrange themselves to the left of the 
imaginary axis. This phenomenon, coupled with the negative sign 
attached to the oscillating terms in the Riemann formula, will be 
enough to offset the impact of the first three negative terms 


5 mer (x?) = : “lit (a) = : Ti (xt?) 


in the Riemann formula, thereby resulting in a(x) > Li(«); ie., the 
prime number theorem undercounts the number of primes. This be- 
havior of the complex zeta zeros is remarkably subtle. 
4.7. Origin of the Riemann explicit formula 
We start with the Euler product. We have 
log ¢(s) = — S/ log (1 — p~*) 
Pp 
1 —ks 
= 277 


p k>l 


lI 
oom 
8 
8 
i 
Q 
YQ 
S 


s(MJ) (s) 
so that 
log ¢(s 
(Ma)(s) = CESS) 
where M is the Mellin transforml| 


‘At the time when Riemann was using the Mellin transform and its inverse, the 
Finnish mathematician Mellin was five years old! Clearly a precocious child. 
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We apply the inverse Mellin transform with c > 1: 


(4.11) ie = = _ 286 (9) neds 
If 
s=5 ti, 
then 
(4.12) E(2) = 5a(9— 1)a-*7T(3/2)¢(s) 


is an entire function of z. Its zeros are given by 

S=pPp, z= -i(p—1/2) =r. 
They lie symmetrically about the origin and are real if the Riemann 
hypothesis is true; and 


(4.13) E(z) = €(0) [a - 27/7’). 


T 


From (4.12) and (4.13) we obtain 
1 
(4.14) log ¢(s) = —log(s — 1) — log’ (1+ 8/2) + gologm + log €(0) 


+ SJ log(1 + (s — 1/2)?/r7) 
and Riemann substitutes from (4.14) into (41) and evaluates the 
terms of the resulting series separately. The dominant term of the 
Riemann formula results from the term — log(s — 1); see §2.12]. 
For the explicit evaluations, see Edwards §3.7]. 


It is interesting to compare the origin of the von Mangoldt formula 
with the origin of the Riemann formula. We recall: 


the von Mangoldt explicit formula comes from the inverse Mellin 
transform of the logarithmic derivative of the Hadamard product for 
the zeta function. 


For comparison, the Riemann explicit formula for J(x) comes 
from the inverse Mellin transform of s~' times the logarithm of the 
Hadamard product for ¢(s). This naturally requires term-by-term in- 
tegration. The detailed computations are in §3.7]. 


Chapter 5 


The prime number 
theorem 


5.1. The Riemann-Ramanujan approximation 


Define 


Re) = S7 A) Lig) 


n=1 
= ie eae ie ae 
= Li(z) Li(a’/*) Li(a’/”) Li(a’/°) 
2 3 5 
1 1 1 
#3 ite) — ; -Li(a/7) + TT Lite ™ 
1 1 
oie Li(a/14) Ts Li(a!/43) reuse 


If, in the Riemann explicit formula, you ignore all the oscillating 
terms, then you obtain the approximation 


which Riemann, in his memoir [R], compares favorably with the 
(Gaussian) approximation 


w(x) & Li(a). 
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Ramanujan’s assertion p. 26] is that 
n(x) — R(x) = O(1), 
ie., that m(a) — R(a) is bounded. 
According to Hardy p. 26], this formula is false. Hardy is full 


of admiration for Ramanujan’s work; nevertheless, this was a famous 
occasion when Ramanujan’s intuition let him down. 


Overall, R(x) is a good approximation to (x). However, it does 
not track the jumps of (x). Furthermore, when a gap occurs in the 
primes (as, for example, between 23 and 29), R(a) keeps increasing 
even though (a) remains constant between primes. 


If we add to R(x) the first few oscillatory terms in the Riemann 
formula (4.10), then we get a new function that very closely matches 
the jumps and irregularities of a(x). When «x is prime, this new 
function increases by about 1 near «, so, in effect, it “knows” where 
the primes are. And when a gap occurs, this new function tends to 
level out, again emulating the behavior of (a). This phenomenon, 
whereby the Riemann formula hugs the function 7(x), has to be seen 
to be believed. 


To see this for yourself, go to www. demonstrations.wolfram.com 
and enter How the Zeros of the Zeta Function Predict the Distribution 
of the Primes in the search tab. If you download the Wolfram Player 
app, then you can vary the parameters. Even with the first 100 pairs 
of zeta zeros, this phenomenon is clearly visible. Try it for yourself! 


5.2. Proof of the prime number theorem 


The prime number theorem is an asymptotic statement. Recall that 
f(x) ~ g(a) as x > co means that 


lim f(a) 


too g(x) 
Theorem 5.1 (Prime number theorem). We have 
m(x)~ Li(w) as «00. 


This theorem was proved (independently) by Hadamard and 
de la Vallée Poussin in 1896. 
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We will give Newman’s proof of the prime number theorem. This 
proof is short and efficient. We will base our exposition on that of 
Zagier [Zag]. Zagier’s exposition is concise and elegant: here we will 
include more details. In order to make this section self-contained, we 
will review and recall certain results from earlier sections. 


We present the argument is a series of steps. Specifically, we 
prove a sequence of properties of the three functions 


C= TS, 06) = OZ 9) = Lose 
n=1 Pp psx 


with sE€C,r@ER. 


Lemma 5.2. We have 


for R(s) > 1. 


Proof. We already proved this; here is a précis of the proof. From 
unique factorization and the absolute convergence of ¢(s) we have 


G(s) = = S) (2™8"8---)-° 


T2,73,-.-20 


for R(s) > 1. 


Lemma 5.3. The function 


extends holomorphically to R(s) > 0. 
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Proof. We already proved this too; here is a précis of the proof. For 
R(s) > 1 we have 


1 1 ae 
_ = Sig as a) 
“@-y = La-f ae 
co n+1 
1 1 
- vf" (4-a)@ 
n nm av 
The series in the last line converges absolutely for R(s) > 0 because 


n+1 n+1 x 
1 aE du 
LO Ge-a)@| = bf [ane 


s 
max a 
aXn<u<n+1 | ust | 


IA 


[s| 
nk(s)+1 


by the mean value theorem. 


Lemma 5.4. We have U(x) < (log 4)z. 


Proof. It is enough to consider the natural numbers. For n € N we 
have 


g2m+1 = (1 op: yess 


2m+1 2m+1 
T+e--+ + +---41. 
m m 


Let M = Ce be the central binomial coefficient. Then we have 


I 


2M < 297+) —» M <2?” —>» log M < 2mlog2. 
Now consider 


0(2m +1) —-V(m+1) = S- log p. 
m+1<p<2m+1 


Each such prime occurs in the numerator but not the denominator of 


w= ae 
m 
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so the product of such primes divides M: 
p|M 
m+1l<p<2m+1 


so we have 


a logp < log M 
m+1<p<2m4+1 
and therefore we have 


0(2m + 1) — V(m+ 1) < 2mlog 2. 


Now use induction, first taking n = 2m +1, an odd number: 
Wn) = VQ@m+1)—V(im+1)4+0(m+1) 
< 22mlog2+ (m+ 1) log4 
= mlog4+(m+1)log4 


n log 4. 
If n is even, then we have 
Vn) = V(n-1) 

< (n—-1)log4 


< nlog4. 


Definition 5.5. The functions f and g are as follows: 


fy: = “Oy (c eR), 
g(s—1): = Ps) — — for (s) > 1. 


Lemma 5.6. We have that g(s—1) is holomorphic for R(s) > 1 and 
g(0) is finite. 


Proof. For R(s) > 1, the convergent product in Lemma [5.2] implies 
that ¢(s) £0 and that 


G(s) log p 


C(s) =a ia 


= @(s)+ S- we i" 
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The final sum converges for #(s) > 1/2, so this and Lemma[5.3]imply 
that ©(s) extends meromorphically to #(s) > 1/2, with poles only at 
s = 1 and at the zeros of ¢(s). However, ¢ is free of zeros on R(s) = 1, 
by Theorem [2.8] So the only possible pole of ®(s) is at s = 1. 


By the elementary Exercise 1.7, the left-hand side of the above 
equation has a simple pole at s = 1 with residue 1. It follows that 
®(s) — 1/(s — 1) has a removable singularity at s = 1. The same is 


true for 5(s) i 
s 

1)= 

OS Tee oa 


so that g(s — 1) is holomorphic for #(s) > 1 and g(0) is finite. 


Definition 5.7. The Laplace transform is defined as 


(LF)(s) = | * e-* F(t)at 
under certain conditions on the ere F and the variable s. 
Lemma 5.8. We have 
(Lf)\(s—1)=g9(s—1) for R(s) > 1. 


Proof. For §t(s) > 1 we have 


B(s) = >> 28? 
7 Pp 
” ir O(a) 
= = 
°° O(x) 
= 7 peri 
= 7 —sty t d 
[ e (e)dt 
and 
ep(s—1) = fee playa 
0 
= —st.9(ef\dt — —(s-ltq 
| e (e)dt | e€ t 
P(s) 1 
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Lemma 5.9. We have (Lf)(0) is finite and equals g(0). In other 
words, the basic Laplace transform equation in Lemma [5.8] remains 
true when s = 1. 


Proof. For T > 0 set 


7 
or(2)= ff er*tat 
0 
This is clearly holomorphic for all z. We must show that 


lim gr(0) = g(0). 


T—00 
Let R be large and let C be the boundary of the region 
{z€C: |z| < R,R(z) > —d} 
where 6 > 0 is small enough (depending on R) so that g(z) is holo- 
morphic in and on C.. Then 


9(0) ~ ar(0) = 5 f (ale) ~ ar(e)je (+ 2/R)S 


Omi 


by Cauchy’s theorem. On the semicircle Ch = CM {R(z) > 0} the 
integral is bounded by 2B/R?, where B = max;50|f(t)|, because 


2) -ortal = | f F(teat 
< B | le **|dt 
ET 
Be-®@T 
OES 
and 
2 
2T a 1 — pR(z)T , 2F(z) 
(0d) eon 


Hence the contribution to g(0) — gr(0) from the integral over Cy 
is bounded in absolute value by B/R. For the integral over C_ = 
CN{R(z) < 0} we look at g(z) and gr(z) separately. Since gr is 
entire, the path of integration for the integral involving gr can be 
replaced by the semicircle C1 = {z € C: |z| = R,R(z) < O}, and 
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the integral over C’_ is then bounded in absolute value by 27B/R by 
exactly the same estimate as before since 


T 

i: f(t)e “dt 

B ge le~**|dt 
—R(z)T 


“RR 
Finally, the remaining integral over C_ tends to 0 as T — oo because 
the integrand is the product of the function g(z)(1 + 27/R?)/z, which 
is independent of T, and the function e*7, which goes to 0 rapidly 
and uniformly on compact sets as T — oo in the half-plane R(z) < 0. 
Hence 


lor (z)| 


IA 


(R(z) < 0). 


lim sup |g(0) — gr(0)| < 2B/R. 
T-00 


Since R is arbitrary, this proves the theorem. 


Lemma 5.10. The following is a convergent integral: 


[ a c 


Proof. We note that the function f is bounded by 1, for 


t 
py = 29 1 <toga-1 <1 
€ 


by Lemma[5.4] Furthermore, we have, by Lemma[.9] 
g(0) )(0) 


= (Lf 

= | toa 

ag 7 
0 

-[° -—2 


Lemma 5.11. We have 
Our) ~ a. 
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Proof. Assume that for some \ > 1 there are arbitrarily large x with 
0(x) > Ax. Since V is nondecreasing, we have 


AL —t AX _ 
/ Ht) dt > if AL t 
x ? x i? 


Xr 
—t 
- 
i 


> 0 


for such x, contradicting lemma|b.10} Similarly, the inequality J(x) < 
Ax with A < 1 would imply 


O)-ty © a —t 


< dt 
Ax ? AL ? 


< 0, 


again a contradiction for \ fixed and x big enough. 


To finish the proof of the prime number theorem, we note that 


Ux) = S- log p 


psu 
< S- log x 
psu 


= (x) logz. 
Rearranging, we obtain 
m(x) log x 
5.1 1< ————.. 
oa = 9G) 


Given 0 < 6 < 1, we have 


U(x) > S- log p 


al~S<p<a 


ye (1 — 6) loga 


= (1—6)(a(@) — m(a*~*)) log x 
> (1 — 6)(a(x) — 21~®) log x. 


IV 
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Rearranging, we obtain 


m(a)logz — logr 4_5 1 
5.2 < ; ee 
oe) ia) 8G)” is 
We have 
See gi-6 = nee +0 as x2 
= : ; oe 
O(x) V(x) =x? 
by Lemma |5.11 
Combining (5.1) and (5.2), we get 
m(x)logx — logr 4_5 1 
1 < . a. 
= oe) ~ ee)” T1=6 


Given € > 0, choose 6 > 0 for which 1/(1— 6) < 1+€/2. Then, there 
exists 29 such that 


m(a) log x 
W(2) 

€ € 

<54+(1+5) 

=1+6 


L>r%o => 1< 


O(a) x 
logx loga 


(a) ~ 
again using Lemma Finally, we have 


ee 
log x Nz) 


by the elementary Exercise 6.1. This concludes the proof of the prime 
number theorem. 
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If y = (a), then we have, successively, 
y~ x/logz, 


It 
ae > 1 
x 
logy + loglogz—logr ~0 as tum, 


as &—> OO, 


logy loglogx 

log x log x 
logy 
log x 


>loas rOM, 


>1loas tom. 


But we have (ylogx)/a — 1; therefore, 


y logy 
x 
as x + oo. Let = pp, the nth prime. Then y = m(p,) = n and so 


> 1 


nlogn 


Pn 
as n — oo. This establishes the following asymptotic statement: 


> 1 


Dn ~ nlogn. 


We may therefore state that the average spacing between primes 
around n is logn. Now multiply n by 10: the average spacing between 
primes around 10n is log 10n = logn + log 10. Therefore: 


Theorem 5.12. If you increase a natural number by a factor of 10, 
then the average spacing between primes will increase by log 10. This 
increase is independent of the chosen number. 


The prime numbers are “irregularly spaced”, but Theorem 5.12] 
pinpoints a certain uniform behavior underlying this irregularity. This 
can be regarded as the down-to-earth interpretation of the prime num- 
ber theorem. Note that log 10 + 2.3. Upon multiplication by 10, the 
average spacing between primes will increase by 2.3. 


Chapter 6 


Oscillation of 7(x) — Li(x) 


6.1. Littlewood’s theorem 


As preparation for the statement of Littlewood’s theorem, we need 
some notation. Let 7(x) denote the number of primes less than or 
equal to x, and let Li(xz) denote the logarithmic integral. The notation 
f(a) = Q4g9(x%) means that 


lim sup f(x)/g(x) > 0, lim int f(x)/g(a) <0. 


There was, in 1914, overwhelming numerical evidence that a(2) < 
Li(x) for all x. In spite of this, Littlewood announced that 


n(x) — Li(x) = 14 (a"/? (log x)~* log log log 2). 


This implies that (x) — Li(x) changes sign infinitely often. 

From a qualitative point of view, the picture which emerges is 
as follows. The average spacing between primes around z is log. 
We recall that we have defined the Skewes number = to be the first 
crossover. As we approach the first crossover =, the primes will start 
to get a little more crowded together, and the average spacing will 
be a little less than log x, sufficiently for w(x) to exceed Li(x) at the 
first crossover. As we move away from the first crossover, the primes 
will gradually thin out, until w(x) is less than Li(x). Then, very 
gradually, the primes will start to get a little more crowded together, 
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sufficiently for (x) to exceed Li(a) at the second crossover. This 
remarkable process will repeat itself infinitely often. 


In the context of the Riemann explicit formula, there will be a 
parallel phenomenon for the zeta zeros. The relevant contribution of 
the complex zeta zeros will be 


-25°R(Li(e*)). 


The issue is one of reinforcement, as opposed to cancellation, 
among the oscillating terms. Around the first crossover =, the terms 
in this infinite series will reinforce one another, in such a way that 
the sum of this infinite series will be sufficiently positive to dominate 
the three large negative terms 


dia a/ay_ Lye ty Ly ats 
5Li(x ) ghi(x ) phi(x ) 


with which the Riemann formula begins. In that case, 7(x) will exceed 
Li(z). This parallel phenomenon for the zeta zeros surely deserves 
further research. 


This chapter contains two oscillation theorems, the first due to 
Littlewood {Li}, the second due to Lehman [Leh]. We will outline a 
proof of Littlewood’s theorem and provide a refinement of Lehman’s 
theorem. 

In Lemma [6.2] below, we have set p = 1/2 + 77 where p is a 
complex zeta zero. 


The method is to average the difference (wu) — wu over a chosen 
interval [e~°x, ez]. If this value is positive, then there must be a 


point within the interval where 7(u) exceeds uw. 


Definition 6.1. This average is denoted as follows: 
1 eax 
F(a,6) := Eres Le —u)du. 


The endpoints of the interval are chosen so that the resulting 
formulas work out nicely. The positive real number 6 is fixed. 
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Lemma 6.2. If the RH (Riemann hypothesis) is true, then we have 


F(x, 6) = 21/2 S- sin(76) sin(y log x) 4 O(ax1/) 
ore Mee Y 


uniformly for x > 4,1/(2x) <6 < 1/2. 


Proof. For this lemma, and the theorem which follows, we will out- 
line the proof. Further details may be found in §15.2 of the book 
by Montgomery and Vaughan and the online class notes of 
Titichetrakun [TT]. 


The starting point is the von Mangoldt explicit formula (£1): 


ue 1 
= ] 1 . 
w(u) =u S- oO os ( =) log 2a 


p 


Integrating from 0 to x, we obtain 


ett 


[ ()- wau = Spry (oem +011) 


for x > 2. We replace x by e*5x and find the difference to obtain 
eo(et+t) _ e—d(e+1) 


26p(p + 1) 


gf Ol). 
We appeal to RH, and rewrite this term as 


6 ei(pt!) — e-d(pt1) 
1/2 74 O(1). 
: (sna) ye 26p(p + 1) ee 


Elementary estimates, as in pp. 477-478], lead to the expression 


_7/2 sinyd 2" O(x'/2 
x yy 6 oa (x*!*). 


By combining the contributions of y and —7, we obtain 


onl? S- sin(76) sin(y log x) 4 O(a'/?) 


ae ar) a 


as required. 
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Theorem 6.3 (Littlewood). As x > oo, 


(6.1) w(x) — & = 04(x/? log log log x), 
and 
(6.2) n(x) — Li(x) = Q4(x!/* (log x) ~! log log log 2). 


Proof. Let N denote a large integer and let 6 = 1/N so that 


F(2,1/N) = —22¥/? S~ are ae me ”) 4 O(@'), 
y>0 


Now the terms 


sin(y/N) sin(ylog x) _ 
{ yi N y ae of 


in this infinite series oscillate in sign, sometimes positive, sometimes 
negative, so it is hard to estimate the sum of this infinite series. 


What if we could find a positive replacement for each term? What 
if many of the terms had modulus close to 1? Then we would have 
a chance of estimating the sum of this series. A theorem of Dirichlet 
Lemma 15.10] allows us make such a replacement. 


Let N be a large integer and set T = Nlog N. Consider those 
numbers y(log N)/27 for which 0 < 7 < T. Define 


K := N(T) X TlogT. 
By Dirichlet’s theorem, there exists an integer 1 <n < N *® such that 


|| = tog | <1/N 
27 


for all0 <7 < T. Take « = N"e*!/" ,§ = 1/N and use the inequality 


| sin 27@ — sin 278| < 27||a — | 
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where ||z|| denotes the distance from x to the nearest integer. We 
obtain 


1 -y/N 
|sin y log x + sin(y/N)| < 2a | 
1 
_s y(nlog N +1/N) + y7/N 
27 
ylog Nn 
=O 
20 
27 
Pena 
~— N 
for0<7<T. Since 
sin(y/N) ! : 
— ft". =| < (log N 
a Y/N ¥ ( ) 


y 


and yop 1/7? «T~*logT <1/N, we infer that 


F(x,1/N) = ¥a/2N- 1S (max) + O(x!/2), 


y>0 


The sum over ¥ is < N log N, so we get 


F(a,1/N) = +221/? log N + O(x'/?). 


But « < NN“ el/N and N(T) x TlogT = N(log N)?, so that 
log log 2 < N(log N)?, 
and hence 
log N > (1 + 0(1)) log log log x. 
The left-hand side in Lemma [6.2] is simply the average of 7(u) — 


over a neighborhood of x. Since « > N and N is arbitrarily large, 
we have (6.1). 


As for we note that if RH holds, then and ( are 
equivalent ee view of Theorem 13.2 in [MV]. If oa is — then 
Theorem 15.2 in gives a stronger result. 


According to Theorem 34], we have the stronger results 


; p(x) - 2 
l > 1/2 
poe Bue log log log x ~ / 
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and 

lim inf ea < 
zoo 21/2 log log log x 

with similar inequalities for 1(a) — Li(z). 


1/2 


6.2. Lehman’s theorem 


Lehman chose to work with the explicit formula for (a) —Li(x) 
as originally conceived in Riemann’s paper. The goal is to develop a 
computational formula for the difference, which depends, to sufficient 
accuracy, on only a finite number of zeta zeros. It is easy to see that 
for any positive function K(y) and any interval (a,b), 


b 
/ K(u){a(y) — Li(y)} dy > 0 


implies that 7(y)—Li(y) is positive within (a, b). Lehman showed that 
with a suitable choice of K(y), the computation of this integral only 
required the first 12,000 zeta zeros to yield a dramatically improved 
upper bound for the Skewes number. Since Lehman’s publication, 
other mathematicians have used his formula with some refinement to 
yield the best upper bound that is currently known. 


Littlewood’s method provided, even in principle, no definite num- 
ber X before which 7(x) — Li(x) changes sign. 


In the course of the 20th century, successive numerical upper 
bounds were found by Skewes conditional on RH, Skewes 
unconditional, Lehman [Leh], and te Riele [tR]. In the course of the 
21st century we have the contributions of Bays-Hudson [BH], Chao- 
Plymen [CP], and Saouter-Demichel [SD]. 

For Littlewood’s own account of the discovery of the Skewes num- 
bers, see pp. 110-112]. 

The smallest value of x with w(x) > Li(x) will be denoted =. It 
is known that 

LO ee 

We now explain the main idea in [Leh]. Lehman’s theorem is an 
integrated version of the Riemann explicit formula. His method was 
to integrate the function u +> 7(e”)—Li(e”) against a Gaussian kernel 
over a carefully chosen interval [w — 7,w +]. The definite integral 
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so obtained is denoted I(w,7). Let p = 1/2 + iy denote a zeta zero 
with y > 0 and let 


iyw 
H(T,w) = 2 SY. “eV Pe, 
0<y<T 

The a in this formula is related to the kernel chosen. Lehman proved 
the equality 

Iw,n) =-1+ H(T,w) +R 
together with an explicit estimate |R| < «. This creates the inequality 

Twn) > H(T,w) - (1+). 
The problem now is to prove that 
(6.3) A(T,w) >1+e. 


If holds, then I(w,7) > 0 and so there exists 2 € [e®~", e? 7] 
for which 7(x) > Li(a). In order to establish (6.3), numerical values 
of the zeta zeros with |y7| < T are required. Each term in H(T,w) is 
a complex number determined by a zeta zero. It is necessary that the 
real parts of these complex numbers, which are spiralling towards 0, 
reinforce each other sufficiently for (6.3) to hold. The only known way 
of establishing this is by numerical computation. When T is large, 
this requires a computer. 


The function H(T,w) is an initial part of the series 


wyw 
H(w) := 2K S> ogee, 
0<y ie 


As Rademacher observed in 1956 [Ra], the Riemann hypothesis 
plus Weyl’s criterion imply that, for each w > 0, the sequence 


{exp(iyw) : ¢(1/2 + ty) = 0,7 > OF 
is equidistributed in the unit circle. So we may expect a fair amount 
of cancellation to take place in the series H(w). This may help to 
explain why it is so difficult to find a number w for which H(T,w) 
exceeds 1. 


We reflect, for a moment, on the Weil explicit formula. This is an 
identity between two distributions p. 339], p. 39]. It is 
well established that certain classical explicit formulas follow from the 
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Weil explicit formula, by picking suitable test functions. For example, 
classical formulas for Dirichlet D-series may be derived in this way; see 
Theorem 3.2, p. 340]. We are led to ask whether the Lehman 
formula can be obtained from the Weil explicit formula by picking a 
suitable test function. 


After this speculation, we will simply state the theorem of Lehman. 


6.2.1. Lehman’s theorem. 


Theorem 6.4 (Lehman [Leh]). Let A be a positive number such that 
B = for all zeros p= B+ iy of C(s) for whichO<y< A. Leta, n, 
and w be positive numbers such that w— 1 > 1 and 


2/A < 2A/a<n<w/2. 
Let 


a =a? 
K(y) = a 2 me, 


w+n 

I(w,) := / K(u—w)ue~"/*{1(e) — li(e“)}du. 
wn 

Then for 2me <T <A 

(6.4) a a Se 


0<|y|<T 


eV /20a4 R 


where R admits an explicit bound with leading term 


3 a A(w a nye @—)/8 


Ww — 1) 


6.2.2. Interpretation of Lehman’s theorem. We will consider 


the term 
el Ww 


gay [20 

p 
which appears in Lehman’s theorem. For those p which satisfy the 
Riemann hypothesis we can obtain this term by substituting y for ¢ 
in the following expression: 


eWP/20(1/2 + it) “te, 
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Now the complex number (1/2 + it)~! tends to 0 as t + oo and 
stays in the sector 


{z : R(z) > 0, S(z) < 0} 
and therefore 
G:th et [2814 Jo +4t)-tet 


is a spiral curve, tending to 0 as t > oo. 


If we replace p = 1/2+7¥ by the conjugate zeta zero p = 1/2—74, 
then the term 


et 


aor [20 
p 


is replaced by its complex conjugate. So the sum of this term and its 


conjugate is equal to 
ome (=) gr ee, 
p 


We note that the complex numbers e’™ are equidistributed on 
the unit circle. Now consider the countable set 


eb 5 1 
—y [2a , gi tela orocg s\ 
{ion ¢($+4) obce. 


Owing to the equidistribution phenomenon, it would seem difficult, a 
priori, for the real parts of these points to be sufficiently negative so 
that 


I(w,n) > 0. 


Littlewood’s theorem guarantees that this phenomenon can in- 
deed happen. 


The main conceptual breakthrough in Lehman’s theorem is to re- 
place the finitely many spirals in the Riemann formula by a single 
spiral. The main practical breakthrough is one of computability: the 
expression I(w,n) is computable. For a large value of T, this will 
require a powerful computer. 


102 6. Oscillation of a(x) — Li(z) 


6.2.3. A refinement of Lehman’s theorem. In the literature on 
the distribution of prime numbers, there is a maze of inequalities for 
1(a). 

We have selected the recent inequalities of Dusart [D] Theorem 
5.1]. These important inequalities are very useful for us. 


Theorem 6.5. For x > 4- 10°, 


(6.5) (a) < = (1+ ae cn 


Ing In?x In’ 


nz 
Forx>1, 
x 1 2 7.99 
: < 14 : 
ee mz) S Ina ( netutetis) 


Since we are dealing with large values of x, we will be using the 
first inequality. The following result is a modest improvement on 
, which is itself an improvement on [CP]. 


Theorem 6.6. Under the hypothesis of Lehman’s theorem and if w— 
1 > 44.22, equation (6.4) still holds if the leading term in R is replaced 
by 


- + log 2+ (w+ nye (w—m/2 


Proof. We recall the definition 
1 1 
J (x) := m(a) + gt (a?) ue gt (ail) fis ise 


and the Riemann formula 


ae d 
. log 2 


I(x) = Li(x) — S > Lif?) + / 
p 
valid for x > 2. We have 


1 1/2 1 1/3 1 1/2 1/3 log x 
= = eee = : 
5M (a y+ gm(z Jers 5m (a y+ gn (a ) lon 


(u? — 1)wlogu 


Now we apply Dusart’s inequality (6.5) to the term $7(«1/), and 
the classic bound m(a) < {225 to the term a(a/3), If > 42-1038, 


— lo 
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then we have 
a ee eee ae 2 8 58.56 
= = sa 14 
57 (a )+ gale Jtri< ar bee ee ines 


logx | «1/8 
log 2 | logx 


and therefore 


See, onl be as gil? 2 8 58.56 
= = bree < 14 
gm )+ gf ) ~ log x log x 


1/3 


+ x 


log 2 


Substituting in the Riemann formula, we have 


l/2 2 8 58.56 
( 2 ( ) ys ( ) log x ( log x log? x — 
2 41/3 
= ee — log 2. 
log2” oe 


uw 
2 
a Pig 9 
pe? e og 2. 
Multiplying through by ue~“/?, we obtain 


(6.7) ue */? (r(e) — Li(e")) > —1— 5" ue“? Li(e™) — F(u) 


p 
where 
2 8 58.56 2Qu 
F(u) := ~¥/68 4 log 2- uel? 
(w) ae a ips og 2 - ue 
Recall that 
Ky) =4f—e ow”, 
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We now integrate (6.7) against the function K between the limits 
w—nandw+n: 


w+n 
/ K(u—w)ue~“/*{n(e") — Li(e“)}du 


> . K(u—w) (-1 — So ue“? Li(e™) — rw) du 
(Sl) =- i. K(u—w)F(u)du 
(92) - [ _ Kiwi 
($3) = i _ K(u—w)ue~"/? (x Lie")) du. 


Our new leading term comes from bounding the expression (S1), 
which is given by 


wen 2 8 58.56 2u 
- K bee al —4/6 4 log 2- ue"/? | du. 
[. (u—w) (= + 7 a 5g og 2- ue du 


In this integral, we exploit the fact that ‘ie K(y)dy = 1, and we 
infer that 
2 8 58.56 
|(S1)| < + 
w-n (w—n)?  (w—7)? 
2(w +) 
log 2 


ee 4. oe (w+ nye m)/2, 


From the point of view of numerical analysis, the crucial improve- 
ment in Theorem [6.6] is in replacing the dominant term 3/(w—1) by 
2/(w — 1). 

Estimating (52) is straightforward. The most laborious part of 
the proof comes up in estimating (53), the sum involving the complex 
zeta zeros, some of which may not lie on the critical line R(s) = 1/2. 
Lehman does not assume any properties of these zeros other than the 
fact that if they exist, the magnitude of their imaginary components 
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exceeds 14, which is approximately the value of the imaginary compo- 
nent of the first zeta zero above the real axis. The rest of the proof is 
a steadfast application of Cauchy’s integral formula to the expression 
in the integral. The details of Lehman’s proof will not be repeated 
here; they can of course be found in his article [Leh]. 


Chapter 7 


The prime number race 


7.1. On the logarithmic density 


Before embarking on computation, we may ask how frequently the 
prime number theorem undercounts the primes. To do this, we need 
a suitable concept of density. We follow p. 174] and p. 
3744]. So, given a subset M C {x € R: & > 2}, the logarithmic 
density of M is defined as usual as 


1 dt 


= lim ——~ : 
X00 log X teMn[2,X] t 


d(M) : 


We are here concerned with the set 
Il := {a > 2: n(x) > Li(x)}. 


Then the logarithmic density 5(II) is a measure of how frequently 
the prime number theorem undercounts the primes. 

Recall the linear independence hypothesis (LI), which asserts that 
the sequence of numbers 7, such tha ¢ (5 + in) = 0 is linearly inde- 
pendent over Q. As usual, the Riemann hypothesis is denoted (RH). 

An important development concerning the “race” between 7(2) 
and Li(x) was the paper of Rubinstein and Sarnak [RS]. They have 
shown p. 175], conditional on (RH) and (LI), that 


6(l) =A 2.6 x 10-7. 
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The conclusion is that, as measured by the logarithmic density, 
there is a definite bias towards 7(x) < Li(x). This may help to explain 
why the Skewes number © is so large. 


Since 7(x) counts primes, it is natural to consider the actual 
primes in the race between 7(x) and Li(xz). What can be said about 
the set of primes p for which z(p) > Li(p)? 

Lichtman, Martin, and Pomerance defined the discrete 


logarithmic density of a set M C Rs2 relative to the prime numbers 
as 


1 1 
§'(M) := jim ——— = 
zoo log log x sep p 


if the limit exists. They also introduced the modified limit 


6*(M) := lim i : logp 
00 OB Tot eu p 


and showed that if the modified limit exists, then it is equal to 6’(M). 


Theorem 7.1 ({[LMP]). Let the set II and the number A = 6(I) 
be defined as above. Assuming RH and LI, the discrete logarithmic 
density of II relative to the primes is 


rah =A, 


This result demands a great deal of reflection. If we may rephrase 
the authors’ own reflections p. 3744], we may say that there 
is no conspiracy for the primes to lie in the set of real numbers IT any 
more or less often than expected. 

As Rubinstein and Sarnak say in their article [RS] p. 175], al- 
though the initial segment in which m(z) loses to Li(«) is extremely 
long, the probability that (x) beats Li(x), although small, is still 
palpable. 


Beating Li(z) is the subject of the next section. 
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7.2. Upper bounds for the Skewes number 


The implementation of Theorem [6.4] consists of first choosing the 
parameters A,a,7,w such that the conditions in Lehman’s theorem, 
Theorem [6.4] are satisfied and then deciding how many actual zeros 
should be entered into the calculation. Then the sum 


S 


o<ini<t 


env (22 


is computed over a wide range of w to find heuristic values of « where 
the sum could be expected to exceed 1. After that, calculations with 
the estimates of the error term R are done to prove that (a) exceeds 
Li(x) around this value. 


Lehman’s formula has the potential to derive a good upper bound 
for the Skewes number because numerical data concerning the zeta 
zeros may be fed into it. 


There are many obstacles in refining Lehman’s result—our imper- 
fect knowledge of the influence of the four parameters and the lack 
of more precise knowledge of the zeta zeros, not to mention our sole 
dependence on Riemann’s formula. 


Where does one go from here? Littlewood’s theorem has illus- 
trated well enough the oscillatory nature of the difference 7(x)—Li(z). 
A possible direction for improvement would be to continue to reduce 
the error term R. There is of course always the option of including 
more zeros in Lehman’s formula, but this will most likely not yield 
a great deal of improvement because the dominant term in the error 
comes from the range of the interval we wish to consider. 


In the 21st century, the first authors to refine Lehman’s theorem 
were Chao-Plymen [CP], swiftly followed by Saouter-Demichel [SD]: 


Year near & N by 
2000 1.39822 x 10346 1x 10° Bays and Hudson 
2010 1.39801 x 10°16 1x 10° Chao and Plymen 


2010 1.3971667 x 1024© 2.2107 Saouter and Demichel 


where N is the number of complex zeta zeros used. 
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We quote from Stoll and Demichel p. 2382]: Our analysis 
conducted with up to 200 billion complex zeros gives computational 
evidence that a crossover may exist near 1.397162914 x 10°16. 

At any rate, a powerful computer is required. Some of the issues 
in handling a large number of zeta zeros are described in and 
will not be repeated here. 

In this context, the following elementary result is useful: if « > 1 
and y > 0, then we have 


"tY dt 
Li(x + y) — Li(z) = i bee 
y 
log x 


With this result, we can state: 


Lemma 7.2. Let 7 > 2 be such that r(x) — Li(x) = M > 0. Then if 
y is a real number such that0 < y <M -loga, we have x(a + y) — 
Li(w + y) > 0. 


Proof. Let y > 0. Since the function 7(x) is increasing, we then 
have 


m(x+y)—Li(x+y) = (m(@+ y) — m(z)) + (r(x) — Li(z)) 
+ (Li(x) — Li(x + y) 
__Yy 
log x 


>M 


This shows that if m(e“) — Li(e”) = M > 0, then a(x) — Li(a) 
will remain positive for another [Mu] consecutive integers. 


After a judicious choice of parameters and after computing the 
first 22 million zeta zeros, Saouter-Demichel Theorem 4.1] ob- 
tained the following result. 


Theorem 7.3. There exists at least one value x in the interval 
[exp 727.9513130, exp 727.9513586| 


for which r(x) > Li(x). Moreover, there are more than 6.09 x 10'°° 
successive integers in the vicinity of exp 727.951335792 where the in- 
equality holds. 
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In the same article [SD], the authors display the Bays-Hudson 
region, the Chao-Plymen region, and their own sharp region. 

Stefanie Zegowitz has obtained a modest improvement on 
Theorem [7.3] with the interval 


[exp 727.951324783, exp 727.951346801]. 


Stefanie’s upper bound provides, at the time of writing , the best 
estimate for the Skewes number: 


19 log 10 < log& < 727.951346801. 
We should also note that, conditional on the Riemann hypothesis, 


Chris Smith, in his MSc thesis, has obtained an improved estimate, 
http: //etheses.whiterose.ac.uk/id/eprint/16409. 


Chapter 8 


Exercises, hints, and 
selected solutions 


8.1. Exercises 


1.1. Compute 7(10), 7(20), 7(30), 7(40), 7(50). 
1.2. Write down a sequence of 333 consecutive nonprimes. 


1.3. (Largest known prime) Estimate the length in kilometers of 
the largest known prime number (see www.mersenne. org). 


1.4. Prove that 


(9 = 


ns 
where d(n) is the number of divisors of n. 


1.5. Prove that 


a(n) 


6(s)¢(s-) => 


ns 


where a(n) denotes the sum of the divisors of n. 
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1.6. (An argument of Euler) What is wrong with the argument 


+00 


1 1/z 
a —0? 
Dd? fe aap 


1.7. (Revision of poles) Let f(s) have a simple pole at s = a. 
Prove that 


ass—a. 


1.8. A computer is programmed to add 1 trillion terms of the 
harmonic series + each second. Since the Big Bang, has enough 
time elapsed for the Nth partial sum to exceed 100? [1 year = 3.156 x 
107 seconds.] 


1.9. Suppose that we used the function li(x) instead of Li(a). In 
that case, what is the Skewes number relative to li(x), i-e., the least 
number for which 7(x) > li(x)? 


1.10. (Life of Riemann) What was Riemann like? When did he 
live? What kind of life did he have? 


2.1. (Theorem of Euler, 1748) Prove that 


fi =1/34/5 1/6 =4/ 741/10 1/11 = 1/134 1/14 4-115 
— 1/17 —1/19 + 1/21 + 1/22 — 1/23 + 1/26 — 1/29 +---=0. 


2.2. The Mobius function p(n) is defined by 
w(1) = 1, wp) = —1, w(p™) = 0 (m> 1), wlab) = wl@)u(b) 
with (a,b) = 1. Prove that for R(s) > 1 we have 


C(s)"* = So u(n)/n’. 
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2.3. (Identity of Ramanujan) Prove that 


(9) _ yn 22 
Ges) = 2 


where w(n) = # distinct prime factors of n. 


2.4. (Revision) Define Laurent series, residue. What is the Lau- 
rent expansion of ¢(s) around s = 1? What is the residue of ¢(s) at 
s=1? 


2.5. (Ramanujan) What manner of man was Ramanujan? Why 
did he come to England? Where did his ideas come from? 


2.6. (Millenium problems, including the Riemann hypothesis) 
Look up www.claymath. org/millenniun. 


3.1. (Euler’s constant ) Let 


ty = 14+1/24+1/3+---+1/(n—1) —logn. 
Prove that t, is an increasing sequence. Write 


n-1 
1 k+1 
r= (- eZ) 


k=1 


and prove that t, is bounded above. Deduce that 
141/2+1/3+1/4+---+1/n—-logn — ¥ 
as n — oo. Deduce also that 
(1+1/2+1/3+---+1/n)/logn > 1 
as N — OO. 


3.2. (Euler product again) Prove that 


G(s) _ yn a(n) 
¢(2s) ns 


where q(n) = 1 if n is square-free and q(n) = 0 otherwise. 
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3.3. Let dy(n) := #{(21,...,2%) € N* : 21-+-2~ = n}. Prove 


that 
g(a) = A) 


n> 


3.4. Look up Euler on www.google.com. What were his main 
contributions to pure mathematics? 


3.5. Prove that 


1 
| _ dx = log 2. 
0 log > 


3.6. Evaluate ¢(—5), ¢(—7), ¢(—9). 
3.7. Find a formula for ¢(—2n — 1). 


4.1. Sketch the step function 
y= 5° A(n), l<a< 30. 


n<x 


4.2. Sketch the step function 
y= >_ An), 1l<a< 30. 


n<ux 


4.3. Sketch carefully the step function 
y=v(2), 1l<a< 30. 
What do you notice about this graph? 


4.4. Show that 
a ee 
pp Pee 
where p = 8 + ty. 


(6 cos(7log x) + ysin(ylogx)) 


4.5. Recall from Chapter 4 the function 


J(a) = n(x) + an(a'?) + an(o¥3) Sos dis 
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which features in the Riemann paper [R]. Show that 


J(x) ~ Li(a) 


as TZ — Oo. 


5.1. (Mellin transform) The Mellin transform is 


yx fH fla)\a*- "de. 


Prove that M(1)(s) =1/s if o > 0. 
5.2. Compute the Mellin transform of x. 


5.3. Prove that 
M(z?)(s) = 1/( =p). 


5.4. Prove that 


co 


M(log(1 — 1/2?))(s) = — $5 1/2n(s + 2n). 


n=1 


5.5. In the formula 


a ier u(r alr) 


Riemann replaced J(x) by Li(a) to e the approximate formula 


oe u(r aT), 


Prove the Ramanujan a. eda 


dr(x aS u(r ar), 
dz aa Dae 


5.6. (Ramanujan identity) Prove that 


(3s) _ d(n’) 
(a) = ne 


where d(n) is as in Exercise 1.4. 
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Ramanujan, 1887 —1920. Indian genius: See The Man Who Knew 
Infinity by Robert Kanigel, Abacus, paperback. See also the novel 
The Indian Clerk by David Leavitt. 


6.1. (L’H6pital’s rule) Prove the asymptotic formula 


x 


ss as log x 


as % > OO. 


6.2. (Ramanujan identity) Prove that 


¢(28) 


where d(n) is as in Exercise 1.4. 


¢°(s) ys: d(n*) 


n> 


6.3. (Ramanujan identity) Prove that 


C4(s) _ yn a(n)? 
Cas) = 2 ne” 


6.4. Sketch the graphs of the first three oscillatory terms in the 
von Mangoldt formula 


We) =e — Soe? /p + O(1). 


6.5. Is the von Mangoldt formula a Fourier expansion? If not, 
why not? 


7.1. Prove that the Fourier series for f(x) = a (—17 < a <n) is 
given by 


1 
and sketch the right-hand side. 


1 1 1 
e=2( Zine Fsinde + zsinae—--) 
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Taking x = 7/2, obtain the Leibnitz formula: 
Ae a A A 


7.2. Prove that, forO<a <1, 
sin(27nz) 
— ~==-2 : 
ae orn 
n>1 


Compare this with the von Mangoldt formula. 


7.3. Define a topology on Z as follows: E is open if EF = 0 or 
FE is a union of arithmetic progressions. Prove that an arithmetic 
progression is open and closed. Define 


A=(J4> 


where A, = {n € Z:n=0 mod p}. Prove that Z— A = {—1,1} is 
not open and A is not closed. Deduce that there are infinitely many 
primes. 


8.1. [Revision] Define Taylor series, Laurent series, isolated sin- 
gularity, isolated essential singularity, pole of order n, simple pole, 
residue, zero of order n, simple zero, holomorphic function in an open 
set, meromorphic function, meromorphic continuation. 


Define the function I'(s) and write down its basic properties. 
Where are the poles of ['(s)? What are the residues at these poles? 
Where are the zeros of I'(s)? 


8.2. The Riemann-Siegel function &(t) is defined by 


1 1 a. a 1 it 1 
S(t) =—- (+ -)a74-?T [i+ — -—+it). 


This function is entire, real valued for real t, and even: =(t) = 


How would you use this function to locate the zeros of ¢(s) on R(s) = 
1? 
5! 
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8.2. Hints and selected solutions 


Notation: O(1) means a bounded function, §t(s) is the real part of s, 
and log x and Inzx both mean the natural log. 


1.1. (10) = #{2,3,5,7} = 4, (20) = (10) + #{11, 13, 17, 19} 
= 8, (30) = 1(20) + #{23, 29} = 10, 1(40) = 1(30) + #{31, 37} 
= 12, ((5) = 1(40) + #{41, 43, 47} = 15. 


1.2. 334! + 2,334! 4+3,...,334! + 334. 


1.3. The largest known prime has 24,862,048 digits, see www. 
mersenne.org. At 1,000 digits/meter, this has length 24, 862 meters 
s 25 kilometers. 


1.4. We have 
1 1 
¢7(s) = Pacecy ae 
m n 
=> 1 
= ae (mn)§ 
ze 
Es I+--41 
= Rc 
ym) 
Ns 
since the number of 1’s is precisely d(N). 
1.5. We have 
1 n 
C(s)¢(s— 1) = eee ie 
3 n 
= et (mn) 
= a(N) 
= aa 


since, for a fixed N, the numbers in the numerator are precisely the 
divisors of N. 
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1.6. The series }> z” is made of two series as follows: 


oo jae, 
1-z 


= 1/z 
1/2" = 1 1 
yen = He hac, 


n=1 


and there is no z which satisfies both inequalities. 


1.7. We have f(s) = (s—a)~1g(s) with g holomorphic and g(a) 4 
0. Then 


f(s) _ (4), (= @)"ta"(s) = (5 ~ aos) 
Oe ray ne (s—a)=*9(5) 
ton 
=*( YG) 1. 


This tends to —1 as s > a. 


1.8. The number of terms is 
15 x 10° x 3.156 x 107 x 1012 & 4.5 x 1079 & 4.5 x e739 ww e667, 


But 
ae 
S- — & log N 
n=1 tt 


and so N x e!°°, Not enough time has elapsed. 


1.9. Here we have 7(2) = 1,1li(2) = 0. The Skewes number 
relative to li(x) is therefore equal to 2. This would, as one reviewer 
wryly observed, spoil the fun! 


1.10. Go to 


www-history.mcs.st-andrews.ac.uk/Mathematicians/Riemann. 
html. 
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2.1. Note the pattern. Only square-free numbers appear. Each 
prime has —1 attached to it. Products of 2 primes have +1. 


Now the Euler product formula for ¢(s) implies that for R(s) > 1 


Since ¢(s) has a simple pole at s = 1, [¢(s)]~! has a simple zero at 
s=1. So if the above equation were valid for s = 1, it would say 


o= 


as required. For a further discussion of this formula, see [Ed] p. 92]: 
von Mangoldt provided a rigorous proof in 1897. 


2.2. The Euler product formula gives 


¢(s)" = []JQ-p™*) =(1- 2°) - 3") -5%)--. 
Now multiply out this infinite product and note that the pattern 
formation is the same as in the definition of pu. 


2.3. The Euler product gives 
¢7(s) Ul (i ae 


“Ua 


¢(28) 


This expression is 
[[G+p-°)(+p-+p-7°+---) = [J + 2p +2p-° +42p- 9+.) 


Now multiply out this infinite product. Each term in the resulting 
sum is of the form 
OX Ox KD Qe) 


kis, kos kps = s 
Py Po ---Dr N 


since r = w(N). 
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2.4. The answers are all in Chapter 1. 


2.5. There are several excellent sources for the life and work 
of Ramanujan. For a detailed account, go to the famous book of 
G. H. Hardy [H]. For a briefer but very incisive account, see the 
chapter on Ramanujan in the book by Stephen Wolfram [W)]. Finally, 
you may go to 


www-history.mcs.st-andrews.ac.uk/Mathematicians/Ramanujan. 
html. 


2.6. Go to www. claymath. org/millenniun. 


3.1. Sketch the graph of y = 1/a and recall that 


in(n) = f° =. 


Mark the points x = 1,2,3,...,2—1,n. Now compare the area under 
the graph with the area under the obvious rectangles. It is then clear 
that t, is an increasing sequence. We have 


The sequence t, is increasing and bounded above and therefore has 
a limit, denoted y. Note that 1/(m — 1) is a null sequence. Since 
tn > 7, we have t,,/In(n) > 0 as required. 
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3.2. We have 


This infinite product is 
(14+27-*)(14+3-*)(14+57%)---. 


Now multiply out. In the resulting sum, only square-free denomina- 
tors appear, with coefficient 1. This coincides with the definition of 


q(n). 

3.3. We have 
a Se eee ee 
1O=. 


aT 
1+1+---+1 
=D 


Ss Ss 
Ny Ny, 


dp.(N 
=) 


since the number of 1’s is precisely d;,(V). 


3.4. Go to 


www-history.mcs.st-andrews.ac.uk/Mathematicians/Euler.html. 


3.5. Make the substitution «++ 2” — x? in the integral 
1 
| (log(1 — x))~"dax. 
2e—e2 


The significance of this example is that the integral of a transcen- 
dental function can be a period “by accident”; see 81.1]. 


3.6. Use the functional equation. 


3.7. Use the functional equation and the Euler formula for ¢(2n). 
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4.1.4.3. Each graph is a step function with jump discontinuities 
at each prime and prime power, i.e., at 


x = 2,3,4,5,7,8,9, 11, 13, 16, 17, 23, 25, 27, 29. 
The jumps are 
In 2,1n3,1n 2, In 5, In 7, In 2, In 11, In 18, In 2, In 17, In 23, In 5, In 3, In 29. 


Note that, in (4.2), the value at « = 7 is n2+In3+1In2+1n5+1n7; 
n (4.1), the value at « = 7 is n2+1n34+1n2+1n5. 

In (4.3), the value at « = 7 is n2+In3+lm2+In5+441n7, the 
average of the values in (4.1) and(4.2). 


In (4.1) each step is an interval of the form (a, b]. In (4.2) each 
step is of the from [a, b). In (4.3) each step is an open interval of the 
form (a,b). 


If your graph is accurate, you will notice that 


W(x) ~ 2; 


ie., U(x) will look like, from a distance, the straight line y = x. 

Since the slope of w(a) is roughly 7/4, the graph is trying to tell 
you that the average space between primes around p is log p. This is 
the first intimation of the prime number theorem. 


4.4. The left-hand side is 


oe 
ape” — op 2 epina 
p 


lp| 

2 aot nz 
a v2 me (pe?! ) 

2x8 - iy ina 
= sap geRlb — iye™ 


and we also have 


RB — iy)e"™™*® = Bcos(ylnz) + ysin(yInz). 
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4.5. Start with the fact that 
1 1 
J(x) — r(x) = 5 tail?) seeds ~a(al/*), 
n 


a finite series. Since 7(x) < x, we get 


1 1 
J(x) — W(x) < ee aed ati 
1 ap 1 aye 
< 5 + he 
1/2 
x 
ey eh eee 
e135 
But k is the integer part of i 5 80 we get 


J(x) — r(x) = O(x'/? log 2). 


Then we have 


Jia) wey, (log x)? 
es = Fe +0 Wea ‘; 


Now let « — oo and use the prime number theorem. 


5.1. Let a = 0 in the solution to Exercise 5.3. 


5.2. Let a = 1 in the solution to Exercise 5.3. We get 


1 
M(c)(s) = —— 
ae 
if R(s) > 1. 
5.3. We have 
M(a2*) = x? s—lda 
1 
_ gas co 
a—s]|, 
eo wl 
~ s—a 


if R(s) > R(a). Now let a = p, a complex zero of ¢(s). If the Riemann 
hypothesis is true, then the condition is R(s) > 1/2. 
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5.4. We have 


l| 


M log(1 — 1/2?)(s) 


I 


n 
1 —2n 
=p ate i) 
Sf By 
i n(s + 2n) 
by the solution to Exercise 5.3 with a = —2n, if R(s) > —2. 


5.5. All you have to do is apply the chain rule to Li(x!/"). This 
gives 


es, 17 ey. 1 1/r-1 
Guile )= agi (1/r)a 
a “ er 
xelnax 


The graph of y = m(a) is a step function, with a jump discontinuity 
at each prime p. If you look at this graph from a great distance, you 
will see a smooth curve, and the Ramanujan formula tells you the 
slope of this curve. 


5.6. See the solution to Exercsie 6.3. 


6.1. As x > oo we have 


-1 
lim Liz) = lim ee) 
x/Inz (nie) —* =e ne) 4 
ah 1 
ae (Ina)-! 
=1. 


6.2. We have the identity 
1-2? 14+2 
(i—a (1-2) 
= (1+ 2)(1+ 22 +327 + 427 +---) 
=1+3x+5x? +72? +--+. 
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It follows that 


Cl) +a =e) 
¢(2s) — II bag 
_y-p*) 
= I (1 = ps) 
=[[]>[@k+pr™. 


Note that d(p?*) = #{1,p,p’,...,p?*} = 2k + 1. Note also that d is 
multiplicative; i.e., d(ab) = d(a)d(b) when (a,b) = 1. The last term 


15 


2k1) 2k 
I~ d(p?*)p—** = » ied 7 gs ) 


Pp Pr 


d(n? 


6.3. We have the identity 
1—«? _ lte 
(er 
= (1+30+507 +723 +---)\(1t+a+a7+---) 


= (k+1)?a*. 


It follows that 


4g t= —s\—A 
can Mapa 
oo) 


pps 
=|] (1 —p-s)4 
=|] > eee 


Note that d(p*) = #{1,p,p?,...,p*} =k +1. The last term is 


: d(pt?)? ---d(pkr)? 
ky2 ks __ 1 r 
WXaohies = ee 
=. % d(n)? 
= rie 
Examples of this kind may be found in Chapter I of the well- 
known book by Titchmarsh [Ti. 
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6.4. If the Riemann hypothesis is true, then the envelope of the 
graph for the complex zero p is y = 2./z/|p|. This curve is the 
parabola 


4x 
y= 2° 
lp| 


As we move up the critical line R(s) = 1/2, these parabolas get closer 
and closer to the z-axis, and the oscillations get faster and faster. 


6.5. The following formulation of the von Mangoldt formula, with 
t = log, has already been proved: 


W(x) -— a ie ~~ sin(Ynt) 
a + O(1) = ae ; 


The fundamental frequency is 71 = 14.184725.... The higher fre- 
quencies 72,73,--. are not integer multiples of y; and so this is not 
a Fourier series, but we do have the same idea as in Fourier series, 
namely building up a function by superposition of waves. 


7.1. We seek the Fourier series of the function f(a), of period 27, 
for which f(x) = x in (—7,7). Now f(a) is odd; therefore, A, = 0. 
We have 


1 Tv 
By, = -| xsinnadx 


T HT 
1 x T 1 

= — |—— cosnz + — cos nxdx 
Tw n —r WJ_y 

= ——cosnt 


so we have 


2 
By, = (-1)"tt=. 


For the graphs, see 


http://en.wikipedia.org/wiki/Fourier_series. 
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7.2. We seek the Fourier series of the function g(x), of period 1, 
for which g(x) = x in (0,1), the Sawtooth wave. We have 


1 
Ao =i xdx = 1/2 
0 
and, when n # 0, 


1 
Ay = 2 | x cos(2rnx)dx 
0 


$ 1 I Laxe 
ay [= ad oF sin(27nz) de 
0 0 


27m 27n 


= 0, 
1 
By, = 2 | x sin(2anax)dx 
0 


_ [een ; 7 of cos(2mnx) | 
0 0 


27m 2mm 


and so we have 


sin(27nz) 
z=1/2 >» aa 


7.3. Each nonempty open set is a union of arithmetic progressions 
and is therefore countable. Therefore {—1,1}, which is the comple- 
ment of A, is not open. So A is not closed. Now each A, is the 
complement of the union of finitely many arithmetic progression, and 
is therefore closed. If there are finitely many primes, then A = J Ap 
is closed, contradiction. 


8.1. The gamma function I'(s) is defined by the equation 


e718 


where 
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I(s) has simple poles at k = 0,—1, —2,—-3,... with residues 


(21 
kl 


It satisfies the relation 
I(s +1) 
Ss 


I'(s) = 


and ['(n) = (n—1)!. The gamma function is zero-free; i.e., it never 
takes the value 0. 


8.2. In the formula for =(t) the terms apart from the zeta function 
are zero-free. So S(t) = 0 iff ¢(1/2 + it) = 0. Now S(t) is real valued 
and continuous, so the zeros on §t(s) = 1/2 can be located by finding 
where &(t) crosses the axis, ie., where =(t) changes sign. This is, in 
fact, how the complex zeros of ¢(s) are found. 
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